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FINAL  REPORT 
OEA  Project  No*  2070 

CONCEPT  AND  FEASIBILITY  STUDIES 
OF  MUZZLE  BRAKE  BLAST  SUPPRESSION  DEVICES 
FOR  105  MM  AND  155MM  HOWITZERS 

1 . 0  INTR  OP UCTION 

Muzzle  brakes  are  included  in  the  new  lightweight  105mm 
and  155mm  Howitzers  in  order  to  reduce  overall  weapon  weight* 
The  weight  savings,  in  the  recoil  system,  results  from 
diverting  a  portion  of  the  propellant  gases  that  normally 
discharge  through  the  muzzle,  after  shot  ejection,  in  a  non- 
axial  direction.  Thus,  the  recoil  system  operates  against 
a  lesser  thrust,  resulting  from  the  Impulse  of  the  discharging 
gases,  and  offers  a  reduction  in  the  weight  of  the  recoil 
system. 

However,  as  the  performance  requirements  of  the  brake 
increase;more  propellant  gas  is  diverted  from  passing  through 
the  weapon  in  the  axial  direction.  The  discharge  of  the 
diverted  high  pressure  propellant  gases  produces  bla s t  waves 
which  in  turn  increase  the  overpressure  in  the  crew  area  of 
the  weapon*  At  present.,  it  is  understood  that  4  psi  overpres¬ 
sure  is  the  maximum  level  which  the  crew  can  tolerate  without 
resorting  to  special  protective  devices. 

This  final  report  contains  a  summary  of  the  work  conducted 
by  Ordnance  Engineering  Associates,  Inc.  under  Contract  Nos. 
DA-11-070-AMC-11  and  01-11-070-003-04909  relating  to  muzzle 
brake  blast  suppression  in  105mm  and  155mm  Howitzers.  The 
main  objective  of  the  program  is  attenuation  of  the  blast  over¬ 
pressure  in  the  gun  crew  servicing  area  without  retarding  the 
effectiveness  of  the  recoil  system. 
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During  the  course  of  this  program  the  results  of  the 
following  studies  were  presented  and  when  practical  applica¬ 
tion  was  made  to  the  105mm  and  155mm  Howitzers; 

1.  Analytical  determination  of  overpressure  maps  resulting 
from  firing  Howitzers  with  muzzle  blast  suppression 
devices* 

2.  Analyses  for  evaluating  the  effects  of  design  parameters 
comprising  the  device. 

3.  Techniques  for  computing  the  dimensionless  pressure  and 
mass  rate  of  discharge  -  dimensionless  time  histories 
in  the  suppression  device. 

4.  The  effect  of  time  which  the  projectile  remains  in  the 
suppression  device. 

5.  Design  parameters  of  the  device  which  create  attenuated 
overpressures  in  the  crew  area. 

6.  A  proposed  design  for  the  105mm  Howitzer  blast  suppression 
device  and  the  stress  analysis  of  the  design. 

In  addition  to  the  above,  an  analytical  study  relating  to 
the  directional  effects  of  blast  waves  was  conducted.  A  summary 
of  this  work  is  presented  in  this  report. 


- <oTa> - - - 
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2.0  DETERMINATION  OF  OVERPRESSURE  RESULTING  FROM 


FIRING  WEAPONS  WITH  MUZZLE  BRAKES 

Overpressures  in  the  crew  area  increase  as  a  result  of 
firing  Howitzers  equipped  with  muzzle  brakes.  In  general, 
the  more  efficient  the  device  the  greater  the  quantity  of 
gas  is  discharged  to  the  rear  area;  hence,  increasing  the 
magnitude  of  the  overpressure.  As  an  upper  limit,  an  idealized 
brake  having  a  maximum  efficiency  would  be  where  all  of  the 
propellant  gas  is  discharged  opposite  to  the  line  of  fire. 

This  ideal  system  would  achieve  the  maximum  reduction  in  the 
recoil  force,  but  would  also  produce  the  greatest  overpressure 
in  the  rear  area  of  the  weapon.  The  magnitude  of  overpressure 
in  the  rear  area  under  this  condition  would  correspond  to  that 
encountered  forward  of  the  muzzle  when  the  weapon  is  fired 
without  a  muzzle  brake. 

In  this  section  of  the  report,  a  method  is  presented 
which  enables  one  to  predict  the  overpressure  resulting  from 
the  blast  of  a  muzzle  brake.  In  addition^the  effect  of  the 
performance  index  and  other  brake  parameters  on  the  overpres- 
8 ure  are  given. 

It  is  necessary  for  a  muzzle  brake  to  have  as  a  minimum 
two  ports  in  addition  to  the  shot  ejection  port  in  order  to 
achieve  non-axial  force  balance.  Such  a  system  is  shown  below: 


—  —  — - — - <QE  A>  -  -  —  -  -  - 
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where 


=  ejection  port  area 

=  recoil  compensating  port  area 

0  =  angle  of  normal  to  A^  with  axis 

R  =  distance  from  brake  to  specific  reference  point 

a  =  angle  of  reference  point  vector  with  respect 
to  axis 


2 . 1  Method  of  Computation 


The  overpressure,  Ap ,  induced  by  energy  flow,  E  ,  from 

1/  e 

a  sonic  orifice,  has  been  shown"  to  be  given  by: 


Ap  = 

poro  'To  ‘  1} 

3 

1/2 

<$Ee 

1/2 

=  K  1/2 

3>Ee 

1/2 

144  R  TT  (Y  +  1  > 
o 

\J 

(2.1) 

ambient  pressure,  psi 

specific  heat  ratio  of  ambient  medium 
distance  from  nozzle  exit,  ft. 

energy  of  gas  discharged  from  orifice,  ft. -lb 

directional  function  dependent  on  angle  of  position 
vector,  R,  with  orifice  axis. 


where : 

P, 


R 

E 

e 


If  the  blast  results  from  a  point  source  detonation,  the 

resulting  shock  wave  is  spherical  as  described  by  Eq.  (2.1), 

and  (J)  =  1.0  for  a  given  energy  release,  E  .  Because  the  energy 

e 

flow  from  an  orifice  is  not  instantaneous,  and  is  also  directed 
in  a  conical  fashion,  the  function  (£)  becomes  strongly  depen¬ 
dent  on  the  angle  from  the  nozzle  axis,  as  well  as  a  factor 

which  may  reduce  the  effective  energy.  This  function  has  been 

2/ 

determined  experimentally^  for  the  flow  from  a  recoilless 

Summary  Report.  Investigation  of  the  Effect  of  Blast  from 
Recoilless  Rifles,  Contract  DA-11 -022-0RD-1 22 7 ,  April  1053 
by  B.  E.  Paul  and  A.  D.  Kafadar. 

— /  Ibid. 
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rifle  orifice,  and  may  be  accurately  represented  by: 


(2,2) 


o 


where : 


6 . 25 


0 


83° 


o 


n  —  3*23 

Thi?  can  be  shown  to  account  for  about  60%  of  the  energy 
as  contributing  to  the  shock  structure  developed.  Of  the 
remaining  40%,  some  of  the  energy  contributes  to  the  formation 


o 

of  very  weak  shock  waves  at  angles  greater  than  83  from  the 


orifice  axis.  The  remainder  simply  does  not  enter  the  primary 
shock  formation  due  to  its  time  lag  in  the  flow  phenomenon. 

However,  Eq.  (2.1)  does  not  yet  lend  itself  to  applica¬ 
tion  to  the  muzzle  brake  situation,  since  there  are  at  least 
two  orifices  to  be  considered  in  the  blast  pressure  developed. 
It  is  possible  to  extend  the  argument  to  include  a  multiple 
orifice  system.  In  general,  energy  is  a  scalar  additive 
property.  Based  on  this  characteristic,  and  assuming  R  large 
with  respect  to  the  maximum  or  if ice-to-or if ice  distance,  the 
energy  function  in  Eq.  (2.1)  may  be  replaced  by  the  summation 
of  directional  energies,  (J)  E^,  contributed  by  each  port  flow. 

This  further  assumes  the  flow  to  occur  in  the  same  time  period, 
which  is  of  course  the  case  experienced  with  the  muzzle  brake 
blast  suppression  device. 

Although  the  energy  term,  E^,  is  equal  for  each  set  of 
symmetric  brake  ports,  the  directional  functions,  (J)^ ,  are 
equal  only  for  a  point,  R,  on  a  plane  containing  the  gun  axis 
and  perpendicular  to  the  plane  containing  the  gun  and  brake 
port  axis.  For  the  case  depicted  in  the  sketch,  the 
functions  would  be 
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1  T 


3>o 


1 


A 


e 

o 


n 


0 


1 


e 


a 


> 


e 


2 


e  + 


(2.3) 


whereas  the  energy  flow  terms,  E  ,  would  be  equal,  E  =  E  . 

1  1 

The  effect  of  the  angle  a  may  be  readily  shown,  however,  by 
referencing  the  overpressure  to  the  gun  axis  for  which  a  =  0. 
The  dual  port  condition  may  then  be  defined  by 


Ap  = 


Ap0  = 


(Apq)  G 

0  n  0  n 

(1  -  e  >  ♦  (1  ' 

_ o _ o 

2  a  -  {->■ 

o 


(K00o)1/2  U  (1  -  *->”  E 


V (2.4  ) 


-il/2 


0  1 
o 


where 


energy  flow  from  each  brake  port 


The  angles,  (0^  0^)  are  defined  by  Eq.  (2*3),  which 

makes  G  a  function  of  the  brake  port  angle,  0,  and  position 

angle,  a„  The  function,  G,  enables  an  axial  pressure  study 

Apo  to  be  extended  to  any  position  angle,  a*  Equations  (2.3) 

and  (2.4)  therefore  describe  the  pressure  distribution  in  the 

crew  area  for  a  double  ported  brake.  Additional  ports  at 

o 

angles  0>a  +  83  make  negligibly  small  contributions.  In 
any  case,  the  technique  is  readily  extended  to  include  any 
number  of  brake  ports,  however,  the  energy  flows  will^in  general , 
be  unequal  and  will  then  enter  as  ratios  in  the  G  function. 
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2  •  2  Brake  Index  ConAtraint 


The  primary  purpose  of  the  muzzle  brake  is  of  course  to 


afford  a  particular  degree  of  recoil  momentum  compensation, 
defined  by  the  brake  index,  i,  as 

Ih 

i  =  1  -  —  (2.5) 

g 


where : 


I 


b 


9 


I 

g 


port  shot-ejection  recoil  momentum  contribu¬ 
tion  of  propellant  gases  with  and  without 
brake,  respectively, 


and  the  resulting  total  recoil  momentum  of  the  gun,  1^,  is 


I  =  mV-  <1-1)1  (2.6) 

r  m  g 

where  : 

m  =  projectile  mass 

V  =  muzzle  velocity 

m 

It  can  be  shown  that  the  momentum  contribution  of  the 
flow  through  an  orifice  is  proportional  to  the  energy  flow,  or 


where : 


=  B  C 


f  i 


B  =  proportionality  constant 

=  thrust  coefficient  of  orifice 
E  =  energy  flow  through  orifice 
1^  =  orifice  axial  momentum 

Without  the  brake,  this  is 


(2.7) 


I  =  B  C  E  (2.8  ) 

g  fg  g 


C 


fg 


1.0,  for  flow  from  a  straight  tube. 
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With  a  muzzle  brake,  the  total  energy  flow  out  of  the 
system  is  less  than  E^  due  to  losses  in  the  brake*  Taking 
this  as  a  function,  b,  of  E^,  yields  the  energy  flow  balance 
as 

E  +  2  Ve  =  (1  -  b)  E  (2.0) 

o  ^  i  g 

where : 

E  =  energy  flow  out  of  ejection  port 

t  h 

E^  =  energy  flow  out  one  of  the  i 
brake  port  orifice  pairs 


During  the  period  of  projectile  base  travel  from  the 

muzzle  to  the  brake  ejection  port,  there  will  be  flow  through 

the  areas  A  ,  however,  there  is  essentially  no  flow  through 
1  1  / 

the  projectile  port,  Aq  Its  contribution  to  the  total  flow 

area  is  therefore  reduced,  so  that  the  energy  flow  distribution 
may  be  written  as: 


E 

o 


XAo 

2  *  \Ao 

AK 

2  SA!  +  \A„ 


(l-b)  E 

g 

( 1  -b  )  E 

g 


(2.10) 


Q  =  effective  area  for  flow  system,  <1 
which  is  seen  to  satisfy  the  energy  balance  given  by  Eq.(2.9). 


Taking  the  gun  axial  components  of  momentum  from 
Eq*  (2.7)  the  brake  Impulse,  Ife,  is  given  by 


or 


I.  =  I  -  2  ^1.  cos  0. 

b  o  ^  i  i 

(2.11) 

I  =  BC  E  -  2B  Vc^4  E4  cos  ©4 
b  foo  *~fii  i 


1/ 


Each  set  of  ports  may  further  be  separated  by  partitions. 
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Substituting  Eq.  (2.10)  into  (2.11),  and  taking  the  ratio  of 
the  result  with  I  from  Eq.  (2.8),  yields  the  brake  index 
a  s  : 


i 


cfoX  “  2  2cfl  ai  cos  ei 


where : 


a 


i 


(2.12) 


The  above  equation  represents  the  brake  Index  constraint 
which  must  be  satisfied. 


2 . 3  Comparison  With  Experimental  Data 


Overpressures  computed  for  the  105mm 

Howitzer  XM103  based 

on  Eqs.  (2.4)  and  (2 

; . 10 )  were 

compared  with  experimental  data 

reported^  for  this 

H  owit zer 

when  fired 

with  three  different 

muzzle  brakes. 

Figs.  2.1,  2.2, 

and  2.3 

compare  the 

above  data  and  were 

based  on  the  following  parameters: 

Brake 

5/K 

D  825  9 

F8241 

X 

1 

1 

1 

p 

0 

0 

0 

Eg  ft"lb 

2. 71xl06 

2. 71xl06 

2. 71xl06 

c 

1 

1 

1 

fg 

c  , 

1 

1 

1 

f  i 

ai 

2. 78 

1 . 28 

3.23 

0  ,  degrees 

25 

30 

25 

P0>  PSi 

14. 7 

14.7 

14 . 7 

1 . 4 

1 . 4 

1.4 

<£>o 

6.25 

6.25 

6.25 

,  d e gr ee s 

83 

83 

83 

n 

3.23 

3.23 

3.23 

““  M uzzle  Blast  Measurement  on  Howitzer.  105mm.  XM103E1. 

H.  H  Holland,  Jr.,  U.S.Aray  Human  Engineering  Laboratories, 
APG ,  TM23-62,  Oct.  1962. 

- <OEA> - 
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Fig.  2.1  -  COMPUTED  OVERPRESSURES  COMPARED  WITH  MEASURED 
105MM  HOWITZER  XM103E1  FIRED  W ITH  5/K  BRAKE 
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It  can  be  seen  in  these  figures  that  the  mathematical  model 
used  to  predict  overpressures  provides  sufficiently  good 
agreement  with  the  experimental  data  to  warrant  its  applica¬ 
tion  to  the  evaluation  of  muzzle  brake  blast  suppression 
devices. 


2.4  E_f  feet  of  Thrust  on  Location s  of  Overpi  essure 

The  following  analysis  was  made  to  show  the  effect  of 
expanding  the  discharged  propellant  gas,  and  thereby  obtain 
additional  thrust  to  reduce  the  recoil  forces,  on  the  location 
of  the  blast  overpressure: 


From  Eq.  (2.12)  the  port  area  ratio  corresponding  to  a 
bi-ported  brake  is  obtained  as: 


.  i 


(i  -  1)  +  (1  -  b)  C 


fo 


(1  -  b)  Cfl  cos  6  -  (i  -  1  ) 


(2.13) 


and  from  Eq.  (2.10)  the  energy  passing  through  each  side  port 
Is  given  by : 


- - — (1  -  b)  E 

2  ai  +  A  s 


(2.14) 


The  location  of  a  selected  overpressure,  Ap,  as  obtained  from 
Eq.  (2.4)  is 


R  =  (J> 


1/3 


-ii/3  r 


P  T 
*o  1  o 


1  4  4tt  (yo+l)(Ap)2 


-rl/3 


(1  -  ^)"+  (1-  ^)n 

o  o 


1/3 


Thus  , 


from  Eq. 

dR 

R 


(2.15) 


(2.15) 


- <OC  A> - 
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— - -  'i  r  i  r  >  M', 


from  E  q .  (2.14) 

dE. 


X  da: 
X+  2ai  a] 


and  from  Eq.  (2.13) 


da. 


(1-b)  cos  6 


dC 


(1-b)  C  cos  0  -  (1-1)  ”fi 


Therefore,  from  the  above  obtain 


M  = 

R 


3  <■ 


_L_, 

X+  2ai 


(1-b)  c 03  0 


(1-b)  Cflcos  0  -  (1-1) 


Ac 


f  i 


(2.16) 


which  shows  the  expected  variation  in  location  of  a  particular 
overpressure  as  a  function  of  the  change  in  the  thrust 
coefficient  of  the  recoil  compensating  ports,  Cf±* 

2.4.1  Application  to  105mm  Howitzer 

The  location  of  4  psi  overpressure  for  105mm  Howit¬ 
zer  fired  with  a  biported  muzzle  brake  was  determined  for  the 
following  cases: 


Case  I 


g 


X 

b 

C 


f  i 


2.714xlOb  ft-lb 
1 . 0 
0 
1 

1.2 


Case  II 


2.714x10  ft-lb 
1 . 0 
0 

1.17  5 
1 . 2 


Case  I  corresponds  to  the  condition  of  equal  pro¬ 
pellant  gas  expansion  in  the  axial  and  side  ports  while  Case 
II  represents  an  Increase  in  the  propellant  gas  expansion  in 


ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


14 


the  side  ports  over  that  in  the  axial  ports.  The  results 
obtained  for  these  cases  are  shown  in  Table  2.1  and  Table  2.2 
for  Case  I  and  II,  respectively. 

Comparison  of  the  data  in  these  tables  shows  that 

the  location  of  4  psi  changes  about  3  per  cent  when  the  C^ 

is  increased  from  unity  to  1.175.  Frc>m  Eq.  (2.16)  it  can  be 

found  that  the  expected  change  in  overpressure  location, 

AR 

would  be  about  25  per  cent. 

Fig.  2.4  shows  a  plot  of  4  psl  Isobar  as  a  function 
of  the  port  angle  0  and  the  location  angle  a.  It  is  seen 
that  as  the  port  angle  0  increases,  the  region  of  4  psi 
overpressure  occurs  closer  to  the  muzzle  of  the  Howitzer. 

The  effect  of  the  parameters  \  and  b  were  checked 
by  re-evaluating  Case  I,  above,  for  0.8  and  b  =  0.05. 

Table  2.3  shows  the  results  of  this  study  and  comparison  of 
these  data  with  those  in  Table  2.1  shows  negligible  variation 
in  the  location  of  the  isobar. 

2.4.2  Application  to  155mm  Howitzer 


A  study,  similar  to  that  conducted  for  the  105mm 
Howitzer,  was  repeated  for  the  155mm  Howitzer  for  the  follow¬ 
ing  cases  : 


E 

g 


X 


b 

C 

i 


f  i 


Case  I 

14.5  6x1 06 
1 . 0 
0 
1 

1.2 


Case  II 

14.56xl06 
1 . 0 
0 

1 . 1  75 

1.2 


These  cases  correspond  to  the  condition  of  equal 
expansion  in  the  side  and  axial  ports  (Case  I)  and  to  greater 
expansion  in  side  ports  than  in  axial  ports  (Case  II).  The 
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Table  2.1 


LOCATION  OF  4  PS  I  OVERPRESSURE  FROM  1QSMM  HOWITZER 

FIRED  WITH  MUZZLE  BRAKE _ WITH  NQ.,  ADD  ITI0NAL-E2lP.AKg.lflM. 


Location  of  Overpressure, 

ft 

\G 

Q  X 

0° 

20° 

40° 

60° 

00 

• 

to 

o 

0° 

10.9 

12.  7 

8.  7 

4. 7 

0 

10° 

A 

14. 7 

13.1 

9.2 

5 . 0 

2.4 

20 

12.0 

14.1 

1  0  ♦  5 

7.3 

5.4 

30° 

10.4 

11.0 

13.8 

9. 1 

0.2 

40° 

8. 3 

10.2 

14.0 

11.0 

8 . 1 

50° 

6.3 

8.4 

13.8 

14.0 

10. 2 

ai 

.  75 

.  8108 

1 . 00 

2.00 

CO 

Table  2.2 

LOCATION  OF  4  PS  I 

OVE  RPRESSURE  FROM 

105  MM 

HOWITZER 

FIRED 

WITH  MUZZLE 

BRAKE 

WITH  ADDITIONAL  EXPANSION 

Location  of 

Overpressure , 

ft 

\G 

a  N 

0° 

20° 

40° 

00° 

80.13° 

0° 

10.4 

12. 3 

8.5 

4.0 

0 

10° 

14.3 

12.  7 

9.0 

5 . 5 

2.1 

20° 

12.2 

13.7 

10.2 

7.2 

3. 9 

30° 

10.2 

11.6 

13.5 

0 . 9 

5 . 8 

40° 

8.1 

9. 9 

13.0 

10.8 

7.8 

50° 

0.1 

8 . 2 

13.5 

14.3 

9.8 

ai 

0.01 

0.00 

0. 05 

1 .54 

CO 

- <Q  E  A> - 
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Fig.  2.4  -  LOCATION  OF  4  PS  I  ISOBAR:  105MM  HOWITZER  FIRED 
WITH  MUZZLE  BRAKE  OF  INDEX  1.2 


Table  2.3 


LOCATION  OF  4  PS  I  OVERPRESSURE  FROM  105 MM  H0WIT2EH 

FIRED  WITH  MU22LE  BRAKE  WITH  NO  ADDITIONAL 

EXPANSION  AND  LOSSES  CONSIDERED 


Location  of  Overpressure,  ft 


\e 

a  \ 

0° 

O 

o 

CM 

o 

o 

60° 

78.27° 

0° 

17.1 

12.8 

8 . 7 

4.7 

0 

10° 

14.9 

13.2 

9.3 

5.5 

2. 3 

20° 

12.  7 

14.2 

10.5 

7.3 

5.3 

30° 

10.6 

12.0 

13.9 

9.1 

6.8 

40° 

8.4 

10.3 

14 . 0 

11.0 

8.0 

50° 

6. 3 

8.5 

13.9 

14.5 

10. 0 
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results  are  summarized  in  Tables  2.4  and  2.5  for  Case  I  and 
II,  respectively.  Fig.  2.5  shows  the  plot  of  the  4  psi  isobar 
as  a  function  of  the  port  angle  9  and  direction  of  angle  a. 
These  data  in  Tables  2.4  and  2.5  show  that  the  effect  of 
increasing  the  expansion  as  indicated  above  serves  to  alter 
the  isobar  location  about  3  per  cent.  This  range  of  varia¬ 
tion  is  in  agreement  with  that  computed  for  the  105mm  Howitzer. 

2 • 5  Summary 

The  results  of  the  studies  conducted  on  the  blast  over¬ 
pressure  resulting  from  firing  Howitzers  with  muzzle  brakes 
are  as  follows: 

!•  Overpressure  can  be  computed  with  the  mathematical  model 
presented  with  reasonable  accuracy  (generally  about  20 
per  cent) 

2.  Directing  the  propellant  gases  away  from  the  crew,  with¬ 
out  reducing  the  effectiveness  of  the  brake,  is  a  good 
method  for  reducing  the  overpressures  in  this  region. 

3.  Providing  additional  thrust  from  the  propellant  gases 
as  they  are  discharged  through  the  recoil  compensating 
ports  should  not  seriously  effect  the  magnitude  of  the 
overpressure  in  the  crew  area. 
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Table  2.4 


LOCATION  OF  4  PSI  OVERPRESSURE  FROM  155MM.  HOWITZER 

F.IRED  .Wi.TH  MUZZLE .BRAKE  WITH  NO  ADD  ITI.QNAL  EXPANSION 

Location  of  Overpressure,  ft 


v 

a  \ 

0° 

o 

o 

<M 

o 

o 

60° 

78.27° 

0° 

29.6 

22.2 

15.2 

8.2 

0 

10° 

25 . 8 

22.9 

16.1 

9. 7 

4.2 

20° 

22.0 

24.6 

18.4 

12.8 

9.5 

30° 

18 . 3 

20.8 

24. 1 

16.0 

10.8 

40° 

14.6 

17.8 

24.5 

19.2 

14.2 

50° 

11 . 0 

14. 7 

24. 1 

25.5 

17.8 

LOCATION  OF  4  PSI 

Table  2.5 

OVERPRESSURE  FROM 

1 SSMM  1 

HOWTT7.ER 

FIRED 

WITH  MUZZLE 

BRAKE 

WITH  ADDITIONAL  EXPANSION 

Location  of 

Overpressure , 

ft 

\  0 

_  o 

_  o 

o 

o 

o 

0 

20 

40 

60 

80.13 

0° 

28.8 

21 . 6 

14.8 

8 . 0 

0 

10° 

25 . 0 

22 . 3 

15.8 

9.6 

3.7 

20° 

21 . 4 

24.0 

18.0 

12.5 

6.9 

30° 

17.8 

20.2 

23.5 

15. 7 

10.2 

40° 

14.2 

17.3 

23.9 

18.8 

13.6 

50 

10.7 

14.3 

23.5 

25.0 

17.1 
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3.0  DESIGN  CONSIDERATIONS 

A  study  was  conducted, based  on  the  analysis  given  in 
Appendix  A,  to  evaluate  muzzle  brakes  and  the  effects  of 
incorporating  blast  suppression  provisions.  The  Influence 
of  the  parameters  comprising  the  device  -  internal  volume, 
v^,  port  angle,  0,  total  recoil  compensating  port  area, 

Ap,  pressure,  p,  and  length,  -  are  evaluated  as  to  their 
Influence  on  the  muzzle  brake  blast  suppression  device 
design.  The  factors  which  comprise  the  suppression  device 
design  are  shown  below: 


A 

P 


total  recoil  compensating  port  area 
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3.1  Application  to  105mm  Howitzer 


Figs,  3,1  and  3,2  were  obtained  from  Eq,  (A. 17)  and 
are  based  on  the  following  parameters  pertaining  to  a  105mra 
H  owi tzer : 


0 


b 


0 

g 

Y 

A 


0,5,  ratio  of  average  gas  tempera¬ 
ture  in  brake  to  isochoric  flame 
temperature 

0.6,  ratio  of  average  gas  tempera¬ 
ture  discharged  from  muzzle  to 
isochoric  flame  temperature 

1.24,  ratio  specific  heat  propellant 
gas 

2 

13,4  in.  bore  area 


C  1  =  V  =  Cd’”  =  0.85,  discharge  coeffi¬ 
cients,  muzzle,  brake  ports,  and 
projectile  exit 


C^'  =  1.3,  thrust  coefficient  muzzle 

=  1.5,  thrust  coefficient  brake  port 

’ 1  =  1.1,  thrust  coefficient  projectile 

exit 


F 

K 


V 

m 


323,000  f t-ibAb,  propellant  impetus 
-6  -1 

6.56  x  10  sec  ,  isentropic  flow 
constant 

12  in.,  brake  length 

4,950  psi,  space  mean  pressure  of 
gun  at  shot  ejection 
-  3 

0.645  x  10  sec,  time  for  projectile 
to  traverse  brake 
3 

1,578  in.  ,  effective  total  gun  volume 
1550  fps,  muzzle  velocity 


The  factor  "a”  indicated  as  a  parameter  on  Fig.  3.1  and 
3.2  is  defined  in  Appendix  A  as  follows: 


a 


v 


b 
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It  can  be  seen  in  Fig.  3.1  that  the  port  angle  required 

for  a  performance  index  of  1.2  is  strongly  dependent  on  the 

factor  "a”  in  the  region  A  /A  between  1.25  and  2;  for  A  /A 

P  P 

greater  than  two,  the  dependence  diminishes.  It  is  also  seen 

that  the  dependence  of  the  angle  on  the  A  /A  ratio  diminishes 

P 

for  A  /A  greater  than  three.  Similar  data  for  other  perform- 
P 

ance  indices  are  shown  in  Fig.  3.2.  It  can  be  noted  in  this 
figure  that  for  a  fixed  port  area  ratio  the  port  angle  must 
be  reduced  as  the  performance  index  is  increased. 

The  above  show  that  operation  with  a  large  port  angle 
(Fig.  2.4  and  2.5  show  that  overpressures  in  the  rear  area 
of  the  weapon  decrease  as  the  port  angle  increases  for  a  fixed 
performance  index)  requires  that  the  port  to  bore  area  ratio 
incorporated  in  the  muzzle  brake  blast  suppression  device  be 
large;  i.e.,  for  an  index  of  1.2  A^/A  should  be  approximately 
three  or  greater.  In  these  regions  the  "a"  factor  does  not 
influence  the  port  angle  profoundly,  but  the  maximum  value  of 
"a"  will  provide  for  an  increase  in  port  angle. 

3.1.1  Effect  of  Design  Parameters  on  Overpressure  Field 

The  effect  of  the  design  parameters  on  the  location 
of  the  overpressure  was  evaluated  by  application  of  the  analyses 
presented  in  Section  2.  The  isobars  mapped  in  Figs.  3.3  and 
3.4  correspond  to  the  parameters  indicated  in  Table  3.1  for 
performance  index  of  1.2. 

Table  3.1 

PARAMETERS  USED  TO  MAP  OVERPRESSURE  105MM  HOWITZER 


\An/A 

Brake  Port 

Angle .  degree 

\P 

a  \ 

1 . 5 

3 

Fie. No. 

1/3 

36.0 

61 . 5 

3. 3 

2/3 

39.5 

62.5 

3.4 
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Fig.  3.3  -  COMPUTED  ISOBARS  OBTAINED  AS  A  RESULT  OF  FIRING  105MM 

HOWITZER  WITH  A  MUZZLE  BLAST  SUPPRESSION  DEVICE,  a  = 
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Review  of  the  above  figures  shows  that  the  overpressure  at 

a  fixed  position  is  reduced  for  the  larger  value  of  A  /A 

P 

in  a  region  subtended  by  an  angle  of  approximately  50  degrees 
with  respect  to  the  Howitzer  axis.  Table  3.2  summarizes  the 
data  showing  the  location  of  4  psi.  It  is  apparent  from  the 
results  given  in  this  table  that  the  overpressure  in  the  crew 
area  can  be  reduced , without  impairing  the  performance  of  the 
muzzle  brake  suppression  device,  by  (1)  increasing  the  brake 
port  area  while  maintaining  the  port  angle  consistent  with 
the  specified  brake  performance,  i.e..  Fig.  3.1,  and  (2) 
increasing  the  parameter  "a*" 

3.1.2  Effect  of  Time  that  Projectile  Is  in  Device 

Fig.  3.5  was  prepared  for  the  105ram  Howitzer  based 
on  the  equations  and  factors  considered  in  Section  3.1.  This 
figure  shows  the  port  angle,  as  a  function  of  the  ratio  of 

the  port  to  bore  area,  A^/A,  for  a  fixed  performance  of  the 
suppression  device  and  for  selected  ’’time'*  and  Ma"  factors. 

It  can  be  noted  that  the  port  angle  is  more  sensitive  to  the 
factor  "a”  the  longer  the  projectile  remains  in  the  suppres¬ 
sion  devices. 

It  is  seen  in  Fig.  3.5  that  for  an  A  /A  ratio 

P 

of  3.5,  the  angle  0  will  vary  as  follows: 

Time  Q.  degrees 

(msec )  1  0.3 

a  =  1/3  66.0  63.8 

a  =  2/3  66.9  64.0 

Hence,  when  the  projectile  remains  in  the  device  for  one  msec, 
decreasing  "a"  from  2/3  to  1/3  corresponds  to  1.4  per  cent 
change  in  the  angle  0.  The  corresponding  angle  change  for  a 
0.3  msec  time  interval  is  0.3  per  cent. 
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Table  3.2 


COMPUTED  LOCATION  OF  4  PS  I  OVERPRESSURE.  1Q5MM 

HOWITZER.  SHOWING  EFFECT  OF  DESIGN  PARAMETERS  y 


Angle  with 
respect  to 
Howitzer  axis 

-  (decree) 

Measured 

Locat Ion 
for  Exit 

of 

of 

4  ps  i 
Brake. 

feet 

 .  a  =  1/3 

a 

=  2/3 

A  / A=1 . 5 

D. 

A  /A=3 

A  /  A=1 . 
D 

5  A 

/A  =  3 

10 

9*3 

5.1 

8.8 

4.9 

20 

10.5 

7.8 

10.0 

6.4 

30 

12.2 

8.5 

11.6 

8.1 

40 

12.5 

10.3 

13.1 

10. 1 

50 

11 . 0 

12.1 

11 . 3 

11.8 
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Ratio  port  area  to  bore  area,  A^/A 


Fig.  3.5  -  BRAKE  PORT  ANGLE  VERSUS  Ap/A  FOR  105MM  HOWITZER 
SHOWING  EFFECT  OF  TIME  PARAMETER 
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It  is  seen  that  the  longer  the  projectile  remains 
in  the  device,  the  port  angle  required  for  the  fixed  perform¬ 
ance  index  increases.  (For  a  fixed  muzzel  velocity,  the  time 
that  the  projectile  remains  in  the  device  is  directly  propor¬ 
tional  to  the  length  from  the  muzzle  to  the  exit  of  the 
device.)  Therefore,  to  reduce  the  overpressure  in  the  crew 
area,  the  time  which  the  projectile  remains  in  the  device 
should  be  maximized.  This  results  from  the  observation  that 
the  overpressure  in  the  crew  area  is  reduced  as  the  port 
angle  is  increased,  thus  directing  the  gases  away  from  the 
service  area. 

3.1.3  Variation  of  Mass  Rate  of  Discharge  through 

Ports  of  Muzzle  Brake  Blast  Suppression  Device 

The  analysis  as  considered  in  this  study  on  the 
rate  of  propellant  gas  discharge  is  based  on  one  dimensional 
sonic  flow  through  brake  ports. 

The  mass  rate  of  discharge  from  the  recoil  compen¬ 
sating  ports  can  be  expressed  as  follows: 

”  *  cd"KAP  Vrf  v  (3-1) 


where 


A 

P 


C 

K 


d 


V  I 


p 


m 


T 

o 


y 


total  brake  pgrt  area  (for  recoil  compen¬ 
sation),  in. 

discharge  coefficient  of  brake  port 

isentropic  flow  constant,  sec  1 

space  mean  pressure  in  weapon  at  time  pro¬ 
jectile  traverses  muzzle,  psi 

o 

gas  temperature  in  brake,  K 

isochoric  flame  temperature  of  propellant 
gas,  °K 

dimensionless  brake  pressure 
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Considering  the  pressure  history  in  the  brake  as  developed 
in  Appendix  A,  it  can  be  seen  from  Eq.  (3,1)  that  the  rate 
of  mass  discharge  is  proportional  to  the  brake  pressure; 
other  parameters  remaining  fixed. 


Fig.  3.6  shows  the  dimensionless  pressure  versus 
dimensionless  time  for  the  following  cases  as  obtained  from 
Eqs.  (A. 11)  and  (A. 12). 


Case 


A  /A 
P  .  . 


I 

1/3 

3.0 

1 1 

1/3 

1.5 

III 

2/3 

3.0 

IV 

2/3 

1.5 

The  mass  rates  of  discharge  from  the  brake  ports  corresponding 
to  the  above  cases  are  shown  in  Fig.  3.7. 

It  can  be  seen  in  Fig.  3.7  comparing  Case  I  with 
Case  III  that  increasing  the  parameter  "a”  from  1/3  to  2/3 
will  provide  a  greater  maximum  rate  of  discharge.  However, 
after  the  maximum  is  reached,  it  is  seen  that  the  rate  of 
discharge  is  less  for  the  larger  value  of  "a"  than  for  the 
value  corresponding  to  Case  I.  Comparing  Case  III  with  Case 
IV,  it  is  seen  that  Case  III  having  the  larger  port  area  ratio 
maintains  a  greater  rate  of  discharge  than  does  the  Case  IV 
with  the  smaller  port  area  ratio  (This  persists  for  about 
15  msec  (T  ^  24 ) . 

3 • 2  Application  to  155mm  Howitzer 

Fig.  3.8  was  prepared  for  the  155mra  Howitzer  and  shows 
the  influence  of  the  parameter  "a”  and  the  ratio  of  port 
area  to  bore  area  on  the  port  angle  for  the  indicated  brake 
indices.  The  factors  used  in  establishing  these  curves  are 
as  foil ows ; 


ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


33 


dJtnssaad  a^eaq  es  a  xuo  t  su  aui  f  p 


ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


34 


PRESSURE  FOR  105MM  HOWITZER 


ORDNANCE  ENGINEERING  ASSOCIATES, 


N  C. 


35 


Fig.  3.7  -  EFFECT  OF  DESIGN  PARAMETERS  ON  RATE  OF  GAS 
DISCHARGE  FOR  105mm  HOWITZER 


Y 


0 


b 


0 

g 


1.24,  specific  heat  ratio  of  propel¬ 
lant  gas 

0.5,  ratio  average  gas  temperature  in 
brake  to  isochoric  flame  temperature 
of  propellant  gas 

0.6,  ratio  of  average  gas  temperature 
in  weapon  (after  projectile  passes 
muzzle)  to  isochoric  flame  temperature 
of  propellant  gas 
2 

29.8  in.  ,  bore  area 

0.85,  discharge  coefficient 

323,000  ft-lb/lb,  propellant  impetus 
-  3  —  i 

6.56  x  10  sec  ,  discharge  coefficient 

6,600  psi,  space  mean  pressure  in 
Howitzer  at  time  projectile  clears  muzzle 

-  3 

0.370  x  10  sec,  time  for  projectile 
to  pass  through  suppression  device 

2700  fps,  muzzle  velocity 

3 

6050  in.  ,  total  effective  gun  volume 


It  is  seen  here  for  a  fixed  performance  requirement  that  as 
one  increases  the  port  to  bore  area  ratio,  the  corresponding 
change  in  port  angle  decreases.  The  effect  of  changing  the 
parameter  "a"  from  1/3  to  2/3  is  seen  to  be  insignificant  in 
the  regions  where  the  slope  of  the  curve  is  not  near  infinite. 
It  can  also  be  seen  that  as  the  performance  index  increases 
the  required  angle  decreases.  Comparison  of  these  data  for  the 
155mm  Howitzer  with  those  reported  for  the  105ram  Howitzer 
(see  Fig.  3.2)  it  can  be  seen  that  the  general  shape  of  the 
families  of  curves  agree.  However,  the  following  exceptions 
are  noted: 


1.  The  155mm  data  shows  less  sensitivity  to  varia¬ 
tions  in  the  factor  "a,"  this  is  attributed  to 
the  difference  in  time  which  the  projectile  is 
in  the  brake. 
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Angle  port  axis  with  respect  to  gun  axis 


90 


Fig.  3.8  -  BRAKE  PORT  ANGLE  VERSUS  Ap/A  FOR  155MM  HOWITZER 
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The  155mm  requires  a  slightly  larger  port  area 
ratio  for  a  fixed  port  angle  and  performance 
index  than  does  the  105mm  Howitzer;  this  is 
attributed  to  the  difference  in  the  pressure 
histories  in  the  weapons  after  shot  ejection. 


3.2.1  Effect  of  Design  Parameters  on  Overpressure  Field 

The  isobars  resulting  from  firing  the  155ram  Howitzer 

with  a  muzzle  blast  suppression  device  were  computed  and  are 

shown  plotted  in  Fig.  3,9,  These  curves  are  consistent  with 

the  parameters  considered  for  the  155mm  Howitzers  in  this 

report.  The  overpressures  as  shown  in  Fig.  3.9  are  plotted 

identically  for  both  values  of  the  factor  "a"  considered,  1/3 

and  2/3.  This  occurs  as  a  result  of  the  insensitivity  of  the 

angle  0  to  the  factor  "a"  as  shown  in  Fig.  3.8.  However, 

it  can  be  seen  in  Fig.  3.9  that  in  the  region  subtended  by  an 

angle  approximately  45°  with  respect  to  the  howitzer  axis,  the 

overpressure  is  less  for  the  A  /A  value  of  three  than  for  a 

P 

value  of  1.5. 


These  results  are  in  agreement  with  those  obtained 
for  the  105mm  Howitzer  (see  Figs.  3.3  and  3.4). 

3.2.2  Effect  of  Time  that  Projectile  Is  in  Device 

It  was  noted  in  Section  3.2  that  the  155mm  Howitzer 
shows  less  sensitivity  to  variations  in  the  factor  "a"  than 
did  the  105mm  weapon.  This  was  in  part  attributed  to  the 
difference  in  time  considered,  that  the  projectile  remains 
in  the  suppression  devices.  Fig.  3 . 10 , prepared  for  the  155mm, 
shows  the  port  angle,  0,  as  a  function  of  the  ratio  port  to 
bore  area  for  a  fixed  performance  index  of  the  suppression 
device  and  for  selected  "time"  and  "a"  factors. 

These  data  are  in  agreement  with  those  shown  in  Fig. 
3.5,  that  is;  the  shorter  the  time  the  projectile  is  in  the 
devices  the  less  sensitive  is  the  angle  to  the  factor  "a"  and 

for  the  longer  the  time  the  greater  the  required  angle  for  a 

fixed  A  /A  ratio. 

P 
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3.10  -  BRAKE  PORT  ANGLE  VERSUS  Ap/A  FOR  155MM  HOWITZER 
SHOWING  EFFECT  OF  TIME  PARAMETER 
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3.2.3  Variation  of  Mass  Rate  of  Discharge  Through 

Ports  of  Muzzle  Brake  Blast  Suppression  Device 

Dimensionless  pressure-time  and  mass  rate  of  dis¬ 
charge  dimensionless  time  curves  were  obtained  for  the  155mm 
Howitzer  based  on  Appendix  A  and  Eq.  (3.1).  Fig.  3.11  con¬ 
tains  the  plot  of  the  computed  pressure  curves  while  Fig. 

3.12  shows  the  mass  rate  of  discharge. 

Comparison  of  the  155mra  curves  with  those  of  the 

105mm  (see  Figs.  3.6  and  3.7)  shows  that  the  effect  of  the 

variation  in  the  factors  "a"  and  A  /A  is  identical  for  both 

P 

systems.  However,  it  is  noted  that  the  maximum  values  are 
greater,  which  is  expected,  for  the  155mm  weapon  than  for  the 
105mm  Howitzer. 

3 . 3  Summary 

It  has  been  shown  that  the  overpressure  in  Howitzer  crew 
areas  can  be  attenuated  by  incorporating  the  following  features 
in  the  design  of  the  muzzle  brake  blast  suppression  device: 

1.  Of  prime  concern,  maintain  the  largest  realistic  port 
area  consistent  with  the  required  performance  index 
(i.e.,  for  the  105mm  a  port  to  bore  area  ratio  of  2.5 
would  be  realistic  for  an  0.9  performance  index  while 

a  ratio  of  about  3.5  would  correspond  as  maximum  for  an 
index  of  1.3). 

2.  Of  secondary  importance,  reduce  or  minimize  the  internal 
volume  in  order  to  maximize  the  factor  "a". 

3.  Also  as  a  secondary  consideration,  maintain  a  large 
length  between  the  projectile  exit  and  recoil  compensating 
ports  in  order  to  increase  the  time  the  projectile  remains 
in  the  device. 


ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


41 


o 


2 

O 

tO 

to 

w 

« 

a 

cu 

D 

to 

H 

to 

< 

CQ 

Cfi 

z 

« 

o 

N 

H 

to 

rn 

a: 

* 

0)  w 

O 

0  H 

X 

•H  W 

4->  S 

s 

< 

ss 

0)  0$ 

m 

m  < 

in 

a>  o« 

rH 

rH 

c  55 

o 

0  o 

w 

•H  H 

os 

0)  to 

p 

C  w 

to 

0)  Q 

to 

0 

w 

•H  Pm 

a: 

T3  O 

a 

•*  H 

w 

a 

K  W 

HH 

Pm 

> 

Pm 

w 

W 

1 

Q 

rH 

rH 

• 

CO 

• 

be 

•H 

Pm 

aanssaad  a^aaq  ssa xuo t eu am f p  *  A 


- ; - <OEA> - 

ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


42 


1  2  00. 


o 


- <6  E  A> - 

ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


to 


43 


DISCHARGE  FOR  155mm  HOWITZER 


4.0  MUZZLE  BRAKE  BLAST  SUPPRESSION  DEVICE  -  105mm  HOWITZER 


A  muzzle  brake  blast  suppression  device  was  made  incor¬ 
porating  those  features  which  attenuate  the  overpressure  in 
the  crew  area  for  a  fixed  performance  index  of  1.2.  This 
device  incorporates  essentially  360  degrees  of  discharge  at 
an  angle  62  degrees  with  respect  to  the  weapon  axis.  The 
design  was  based  on  the  following; 

1.  Ratio  of  port  to  bore  area  of  about  three. 

2.  Internal  volume  minimized  consistent  with  an 
overall  length  of  13.5  in.  and  yielding  an  ,,a" 
value  of  0.59. 

3.  62  Degree  port  angle  selected  to  yield  perform¬ 
ance  index  of  1.2. 


4  •  1  Design  Description 

The  design  of  the  360°  muzzle  blast  suppression  device 
shown  in  Fig.  4.1,  consists  basically  of  a  threaded  ring  for 
attachment  to  the  Howitzer  muzzle,  a  conical  shaped  baffle 
with  a  cylindrical  extension  and  a  set  of  four  equally  spaced 
gussets  which  join  the  ring  and  baffle.  The  port  area  is 
located  between  the  baffle  and  ring  in  the  axial  direction 
and  the  gussets  in  the  circumferential  direction.  Design 
values  for  the  device  are  itemized  below: 


a  . 

b. 

c . 
d  . 
e  . 
f  . 


g- 


Howitzer  bore  area  .  . 

Port  area  . 

Entrance  to  port  area 
Port  to  bore  area  ratio 
Internal  volume  .  .  . 

Internal  length  .  .  . 

Parameter  "a"  .... 


A  =  13.4  In.2 

A  =  41.3  In.2 

P  2 
A  =  47.2  in. 

e 

.  A  /A  =  3.08 

P  3 

v.  =  223  in. 
b 

L .  =  9.88  in. 

b 
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h.  Port  angle  . 0  =  62° 

i.  Approx,  weight  of  device  ....  W  =  59  lb 

In  Fig.  4.1  the  device  is  shown  fabricated  as  a  casting. 

As  a  weldment  the  gussets  would  be  modified  to  permit  access 
for  welding  to  the  baffle,  increasing  the  weight  of  the  de¬ 
vice  by  approximately  4  lb. 

The  results  of  the  stress  analysis  indicates  a  maximum 
tensile  working  stress  of  57,100  psi;  with  a  minimum  yield 
stress  of  120,000  psi  for  the  material,  a  factor  of  safety  of 

2.1  is  obtained. 

4 . 2  Results  -  Stress  Analysis  of  360  Degree  Muzzle 

Brake  Blast  Suppression  Device 

The  stress  analysis  of  the  360  degree  device  for  the  105mm 
Howitzer,  involved  the  determination  of  maximum  stress  levels 
and  deflections  when  subjected  to  a  maximum  gas  pressure  of 
5000  psi.  Simplifying  assumptions  were  made  so  that  the  stress 
distribution  at  the  statically  indeterminate  junction  between 
the  gussets  and  baffle  could  be  evaluated.  It  was  assumed  that 

1.  Uniform  stress  distribution  exists  between  gus¬ 
sets  and  baffle  in  the  axial  direction.  As  a 
result  of  this  assumption  the  longitudinal  bending 
moment  at  the  junction  of  gussets  and  baffle  is 
neglected . 

2.  Localized  baffle  deflection  in  the  radial  direc¬ 
tion  at  the  gussets  can  be  approximated  by 
utilizing  the  ring  theory. 

In  addition,  the  principle  of  stress  superposition  was 
incorporated  in  this  analysis  to  find  maximum  combined  stress 
le ve 1 s  . 
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The  method  of  solution  included  evaluation  of  the  max¬ 
imum  radial  disparity  between  gussets  and  baffle  when  separated 
and  gas  pressure  acting  upon  the  baffle  alone.  The  disparity 
is  removed  by  introducing  radial  forces  between  gussets  and 
baffle  and  a  bending  moment  to  the  mounting  ring.  Maximum 
stresses  are  determined  by  methods  of  superposition  of  the 
resulting  stresses.  In  addition,  the  tensile  stress  in  the 
gussets  resulting  from  the  gas  pressure  impinging  upon  the 
baffle  in  the  axial  direction  was  included. 

A  summary  of  the  computed  stress  levels  and  corresponding 
strains  is  given  in  Appendix  B  and  is  presented  in  Table  4.1, for 
the  sections  depicted  in  Fig.  4.2. 

4 . 3  Summary 

Selection  of  the  circumferential  type  port  was  based  on 
establishing  the  conditions  in  the  device  such  that  the  gases 
passing  through  the  port  are  directed  at  the  greatest  angle 
with  respect  to  the  Howitzer  axis  while  maintaining  a  fixed 
performance  index.  This  circumferential  port  arrangement  pro¬ 
vides  the  following  conditions  which  are  consistent  with  the 
above  criteria: 

1.  Port  to  bore  area  ratio  maximized  consistent 
with  minimum  volume  and  realistic  length. 

2.  Internal  volume  minimized. 

3.  Internal  travel  of  the  projectile  while  in 

the  device  maximized  consistent  with  a  reason¬ 
able  length  device  and  by  providing  the  gas 
ports  as  close  to  the  weapon  muzzle  as  possible. 


- <OEA> - 
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Table  4,1 


SUMMARY  OF  MAXIMUM  STRESSES  AND  D ISPLACEMENTS 

FOR  360°  MUZZLE  BLAST  SUPPRESSION  DEVICE 


Loca  1 1  on 

(  See 

Fie.  4.2) 

Maximum 

C  omb i ned 
Stress 
(osi  ) 

Direction 

of 

Stress 

Radial 

D isplaceraent 

<in» .) 

(1  ) 

53 , 000 

Circumferential 

0.0062 

(2  ) 

57,100 

Longitudinal 

0.0005 

(3  ) 

34 , 300 

Longi tud inal 

0. 0005 

(4) 

36 , 200 

Circumferential 

0.0029 

(5) 

16,500 

Longitudinal 

0.0005 

(6) 

19,400 

Longi tud inal 

- 

gus 


le 


Fig.  4.2 


LOCATION  OF  COMPUTED  STRESS  AND  DEFLECTIONS 
FOR  360°  MUZZLE  BRAKE,  105MM  HOWITZER 
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5 . 0  PRESSURE  AND  FLOW  FIELD  EXTERIOR  TO  A  GUN  BLAST 

The  axi-symme trie  acoustical  problem  of  a  gun  blast  has 
been  formulated  as  an  initial  value  problem.  The  exact  solu¬ 
tion  is  obtained  for  a  given  initial  pressure  and  velocity 
distribution.  Special  cases  studied  in  detail  Include: 

(a)  an  initial  velocity  distribution  with  a  semi-infinite 
cylinder,  and  Cb)  a  constant  pressure  distribution  within  a 
finite  cylinder. 

An  asymptotic  representation  for  the  pressure  field  is 
also  developed  to  gain  insight  for  the  far  field.  Several 
plots  for  the  acoustical  pressure  field  are  given. 

The  stydy  of  the  problem  is  presented  in  Appendix  C. 

This  study  was  carried  out  by  Dr.  A.  C.  Eringen  of  Purdue 
University  under  OEA  Purchase  Order  No.  5889.  The  analysis 
presented  was  reviewed  by  OEA  personnel,  especially  by 
Dr.  Nurettin  Y.  Olcer. 
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6  *  0  RECOMMENDATIONS 


The  results  of  these  analytical  studies  have  shown  that 
overpressures  in  the  crew  area  can  be  reduced  without  affect¬ 
ing  the  performance  of  the  muzzle  brake  by  selection  of  those 
design  parameters  which  allow  the  gases  to  be  directed  at 
the  greatest  angle  with  respect  to  the  weapon  axis.  Based 
on  the  results  obtained  under  this  program  it  is  recommended 
that  additional  studies  be  undertaken, 

1.  Verify  experimentally  the  analytical  results 
reported , 

2.  Modify  analysis,  if  necessary,  based  on  ex¬ 
perimental  results  and  prepare  general  design 
manual  for  muzzle  brake  blast  suppression  devices. 
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APPENDIX  A 

INTERIOR  BALLISTICS  OF  MUZZLE  BRAKE 

AND  BLAST  SUPPRESSION  DEVICES 
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APPENDIX  A 


INTERIOR  BALLISTICS  OF  MUZZLE  BRAKE 

AND  BLAST  SUPPRESSION  DEVICES 


A. 1  BRAKE  PRESSURE 


The  pressure-time  history  in  the  muzzle  brake  was  devel¬ 
oped  based  on  the  following: 

1.  Hugoniot  form  describes  the  pressure  history 

in  the  weapon  after  the  projectile  has  traversed 
the  muzzle. 

2.  The  projectile  traverses  the  brake  with  uniform 
vel oc i t y . 

From  the  Hugoniot  Equation  the  space  mean  pressure  in  the 
weapon  after  the  projectile  has  left  the  muzzle  is  given  by: 


p  =  p  (1  +  B  t)  Y_1 
m  l 


(A.  1  ) 


where 


P 

P 


m 


P 

F 


space  mean  pressure 


1/ 


space  mean  pressure  at  the  instant 
the  projectile  is  discharged  from 
the  muzzle. 


(1  +  P) 


fraction  of  heat  loss  in  weapon 
propellant  impetus 


A  =  bore  area 


The  space  mean  pressure  is  the  order  of  90  per  cent  of  the 
breech  pressure  for  the  type  of  weapon  considered  in  this 
study.  However,  to  simplify  the  analysis  and  without 
seriously  affecting  the  results  it  shall  be  assumed  that 
the  ratio  of  space  mean  to  breech  pressure  p/p  is  equal 
to  unity. 
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r 

T 

m 

T 

o 

K 

K 


total  gun  volume 

discharge  coefficient,  muzzle 

ratio  of  specific  heats  of  propellant  gas 

gas  temperature  at  shot  ejection 

isochoric  flame  temperature  of  propellant 
gas 

Isentropic  flow  constant 


X£ 

F 


X±i  1 

2  Y-l 

<-*7> 

Y+l 


1/2 


g  =  gravity  acceleration 


The  pressure  history  in  the  brake  during  the  period  when 
the  projectile  is  present  can  be  obtained  from  the  following: 

Energy  Balance 


c  NT 
P  g 


Pi 


c  NT 
P  g 


+  c  (N-N,  )  T.  + 


:  N.T. 

p  1  b 


(A. 2) 


Equation  of  State 
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-¥*5 


where 


N 

T 

g 


quantity  of  gas  entering  the  brake 
gas  temperature  leaving  barrel 
gas  temperature  in  brake 

isobaric  specific  heat  of  propellant  gas 
isochoric  specific  heat  of  propellant  gas 
discharge  coefficient  brake  ports 
quantity  of  gas  discharged  from  the  brake 
fraction  of  heat  loss  in  brake 
brake  pressure 
bore  area 

length  of  projectile  traverse  in  brake 

t*ime  for  projectile  to  traverse  at 

constant  velocity 

total  port  area  of  brake,  does  not  include 
projectile  exit  area 

interior  volume  of  brake 


C  ombi ni ng 

Eq.  (A. 2) 

(A.  3  ) 

and  (A. 

4  )  one 

can  obtain 

djr 

h 

B1 

( 1  +  bT) 

_2X 

r-i 

O^T^l 

dT  + 

1  +  rT  y  " 

1  +  rT 

t 

(A. 5) 


where 
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A-3 


’J 


f 


f 


9 


0 


2 

1 

g 

b 


a 


b 


1 2FC  .  1  *  KA  t, 

__  jd _ _  p  i 

vb 

1 2FC  , * KA  t, 
a  1 


T 

T 

o 


T 


o 

AL, 


B  t 
11 


During  the  period  after  the  projectile  has  passed  through 
the  brake,  the  basic  equations  become: 

Equation  of  State 

Th 

pbVb  =  12F  (N-Ni  )  t  (A. 6) 

.  o 

and 


The  energy  balance  and  rate  of  gas  discharge  are  as  given  above. 

The  corresponding  differential  equation  resulting  from 
the  combination  of  the  above  system  of  equations  is 

_  _ 

^  +  My  =  ( 1  -a  )  gx  Cl  +  bT  )  Y_1  ,  T2  Tj_  (A. 8) 

where 

m  =  f  <w 
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The  solution  to  Eq.  (A. 5)  is  given  by: 


y  =  C1(l+rT)  - ( 1+rT) 


h 

T 


J 


gl  (1  +  b^  )  ^  1  (l  +  r<£  ) 


h  +  r  -\ 

T 


0£T£1 


(A.  9) 


and  for  Eq.  (A. 8)  by 

-MT  -mT 

y  =  C  x e  -e  g^l 


-a)  IA 


_  2 £_ 

(1  +  b^  )  Y  1  (A. 10) 


T21 

The  Integrals  of  Eq.  (A. 9)  and  (A. 10)  are  obtained  based 
on  the  following  assumptions 


_  g.i_. 

Cl  +  b^  )  Y“1  =  (1-m^),  O^T^l 

m  -  _  b 

Y“1 

This  assumption  Is  justifiable  since  b  Is  small;  the  order 
of  10  ^  and 


(l  +  b^  ) 


?r  .  A 


=  e  °  ,  T>  1 


The  solution  to  the  differential  equations  describing  the  di¬ 
mensionless  pressure-time  history  In  the  brake  are  as  follows: 


®*1 

y(T)  =  -ji 


i+* 

r 


1- (1+rT) 


gl 

m 

f 

Ch  +  r ) 

- 

r 

»  ^ 

( 

h  -i 


0<J£1 


(A. 11  ) 
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-MCT-l  ) 


y  (T)  =  yxe 


+ 


1  -a 

M-n 


-m  T  -M(T-1  )-m 

e  -e 


,  T>  1 

(A. 12) 


a  nd 


=  y(l)  from  Eq.  (A. 11) 


A. 2  BRAKE  PERFORMANCE 


From  the  definition  of  brake  index,  the  following  can 
be  written 


C  O  8  0 

I  +  (i-l)  I 
e ,  e 

(A. 

I 

D 

where 

0  = 

angle  of  gas  deflection,  through 

thr us t 

reversing  ports,  with  respect  to 
weapon  axis 

the 

i  = 

muzzle  brake  index 

I  = 

e 

momentum  produced  by  the  gas  flow 
the  shot  exit  port* 

through 

I  = 

P 

momentum  induced  by  the  gas  flow 
the  brake  reversing  ports 

through 

I  - 

B 

momentum  resulting  from  gas  flow 
muzzle  with  no  brake  present 

through 

The  impulse  as  shown  above  are  obtainable  from  the  following 


I 

e 


I 

P 


r*  co 


p.  dt 


c  f ' A 


J^2Pmtl 


(A. 14) 


Pm1! 


(A. 15 ) 


I 

g 


rc 


Cf,,,  A 


pdt  —  C  "fA  A_p  t. 
c  f  *  iV3  m  1 


(  A  #  1  6  ) 
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where 


Cf  =  thrust  coefficient  corresponding  to 
each  of  the  ports;  muzzle,  brake 
reversing  and  shot  exit 

i 

A1  =  ydT 

O-' 


Substituting  Eqs.  (A. 15)  and  (A. 12)  into  Eq.  (A. 13),  obtain 


cos  0  = 


V 


A 


2  +  (I-l)  cfH* 


A: 


A 

_E 


CA.17) 


f  A 


<  Ai  +  A2> 


The  A  terms  are  evaluated  from  Eqs.  (A.l),  (A. 11)  and  (A. 12) 

as 

£  i 


Ai  =  r 


i  ♦  = 
r 


i  -  * 

1  .  r  -i 


r  ( 1  -  -) 
r 


h  +  r 


1  -  - 

2+r  _  (1+r)  r  -1 

2  r ( 1  -  -) 


(A. 18 ) 
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A. 


1 

M 


(1-a)  g 


*1  + 


m 


1  -m 
-  e 


(A. 19) 


1 

b 


>  ) 
Y+l 


(A. 20) 
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APPENDIX  B 

STRESS  ANALYSIS  OF  MUZZLE  BRAKE  BLAST 

SUPPRESSION  DEVICE  -  105mm  HOWITZER 
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STRESS  ANALYSIS  OF  MUZZLE  BRAKE  BLAST 

SUPPRESSION  DEVICE  -  105mm  HOWITZER 


The  design  of  the  muzzle  brake  blast  suppression  device 
shown  in  Fig.  4.1  was  based  on  the  stress  analysis  given  in 
this  appendix. 

B.l  EQUATIONS  FOR  CONICAL  BAFFLE 

The  circumferential  stress—^  in  a  conical  section  with 
internal  pressure  is  given  by 


o 

PRi 

S1 

t^  cos  y 

where 

P 

= 

internal  gas  pressure, 

ps  1 

Ri 

= 

mean  radius  of  wall,  in 

. 

h 

= 

wall  thickness,  in. 

Y 

= 

wall  slope,  deg. 

The 

corresponding 

radial  displacement  is 

given 

ARX 

= 

si  Ri 

E 

or 

pR  2 

AR1 

1 

Et^  cos  y 

where  E  = 

modulus  of  elasticity,  ps i . 

(B.l  ) 


(  B  .  2  ) 


The  conical  baffle  wall  shape  was  proportioned  so  that 
the  hoop  stress  is  essentially  constant;  the  radial  displace¬ 
ment  under  pressure  loading  is  then  proportional  to  R^.  The 

Roark,  R.  J .  ,  Formulas  for  Stress  Analysis  and  Strain , 

Me G raw  Hill,  N.Y.,  1954,  page  269. 
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gussets  act  to  restrain  this  displacement  giving  rise  to 

inward  forces  at  the  junction  between  gussets  and  baffle. 

For  each  circular  ring  segment  of  the  conical  baffle  the 

following  equations  were  derived  by  superposition  of  180° 

2/ 

loading  of  a  circular  ring-  giving  the  following  equations 
o 

for  90  loading  of  the  ring# 


M  = 

AR1f 


0.136  WRx 


=  -0.060 


(B .  3  ) 

(  B  .  4  ) 


S 


b 


(  B .  5  ) 


where 


M 


M  • 


W 

AR 


« 

1 

« 


maximum  bending  moment  in  circumferen¬ 
tial  direction,  in. -lb 

maximum  unit  bending  moment  in  circum¬ 
ferential  directions,  in. -lb  per  in.  of 
baf  fie  le  ngth 

radial  inward  loading  of  ring  section,  lb 

area  moment  of  inertia  of  ring  section, in. 

radial  inward  deflection  of  ring,  in. 

bending  stress  in  circumferential  direc¬ 
tion,  ps  i 


B.  2  EQUATIONS  FOR  GUSSETS 


The  gussets  were  considered  to  be  short  cantilever  beams 
joined  to  the  mounting  ring  and  loaded  in  the  vertical  plane 
at  the  junction  with  the  baffle.  As  a  conservative  estimate 
it  was  assumed  that  the  mounting  ring  and  gussets  are  not 
affected  by  the  gas  pressure  or  radial  expansion  of  the  howit¬ 
zer  muzzle. 

2/ 

—  Ibid,  page  156,  case  2. 
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Radial  displacement  of  the  gussets  at  the  Junction  with 
the  baffle,  considering  the  unloaded  end  fixed  to  the  mounting 
ring,  is  given  by 


AR 


2 


Wg3  6WI 

3EI2  5A2G 


(B.6) 


where  the  first  terra  represents  the  contribution  from  bending 

and  the  second  terra  represents  the  contribution  from  shear 
3/ 

loading—,  and 

W  =  end  loading  of  gusset  during  radial 
displacement  of  baffle,  lb 

%  -  length  of  gusset  to  mounting  ring, 

neglecting  remainder  of  gusset,  in* 

=  area  moment  of  inertia  of  gusset  at 
support  in  vertical  plane,  in. 

G  =  modulus  of  rigidity,  psi 

Tensile  stress  in  the  gussets  in  the  axial  direction  during 
pressure  loading  of  the  baffle  is  given  by 


S 


2 


PAb 


A 

g 


(B.  7) 


where 


P 

A 


b 


A 

g 


muzzle  brake  pressure,  psi 


net  projected  area  of  conical  baffle 
in  vertical  plane,  in.2 


cross-sectional  area  of  four  gussets 
in  vertical  plane,  in.2 


Bending  stress  in  the  gussets  at  the  Junction  with  the 
mounting  ring  is  expressed  as 


b  2 


(B.  8  ) 


where  Z  =  section  modulus  of  gusset  in  vertical 

pla  ne ,  in. 

3  / 

Ibid,  page  119. 
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B .  3  EQUATIONS  OF  MOUNTING  RING 


The  mounting  ring  is  loaded  at  the  front  edge  with  con¬ 
centrated  bending  moments,  radial  shear  forces  and  longitu¬ 
dinal  tension  forces  located  90°  apart*  Outward  flaring  of 
this  end  of  the  ring  serves  to  reduce  the  radial  load  W  and 
consequently  the  bending  moment  M  in  the  conical  baffle 
section.  As  a  conservative  estimate  it  was  assumed  that  the 
bending  moments  and  shear  forces  transmitted  from  the  gussets 
are  uniformly  distributed  around  the  circumference  of  the  front 
edge.  Then  the  contributions  to  the  radial  displacement  of 
the  gusset  baffle  Junction  by  deformation  of  the  mounting  ring 
is  given  by: 


where 

4> 


KW 

o 


However“{ 


where 


M 

o 


AR3  =  t<f>  *  KWo 


(B.9) 


=  angular  displacement  of  front  edge 
of  mounting  ring,  radians 

=  radial  displacement  of  mounting  ring 
resulting  from  shear  load  W,  in. 


4> 


M  W 

o  o 

2DX2 


(B.10) 


in. -lb/in.  (B.ll) 

^3 


W 

o 


lb/in. 


(B.ll) 


K 


2D\2 


M 

<T  +  ~)  in./lb 
A  Wo 


(  B  .  1 2  ) 


4/ 


Ibid  , 


page  271. 
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X 


per  in 


( B  •  1  3  ) 


D 


a  nd 


3 
V 
t  0 


- 2  T 

2  2 
R  * 

3  3 


E  t . 


12  (1-V  ) 


lb-in. 


=  mean  radius  of  ring,  in. 

=  Poisson#s  ratio 
=  wall  thickness  of  ring,  in. 


(  B  .  1  4  ) 


The  maximum  bending  stress  in  the  front  edge  of  the 
mounting  ring  is  given  by 


S 


b3 


6Ho 


1 . 932  Wo 


lb/in. 2 


(B. 15  ) 


The  tensile  stress  in  the  longitudinal  direction  resulting 
from  longitudinal  forces  transmitted  from  the  gussets  is 
assumed  uniformly  distributed  and  approximated  by 


S 


1 3 


A 


r 


lb/in. 2 


(B. 16  ) 


where  A 

r 


cross-sectional  area  of  ring  in  vertical 
plane . 


B .  4  SUPERPOSITION  OF  DISPLACEMENTS  AND  STRESSES 

B . 4 . 1  Radial  Displacement  at  Gusset-Baffle  Junction 

The  maximum  radial  displacement  of  the  conical 
baffle  without  gussets  attached  occurs  at  the  outer  diameter 
as  given  by  Eq.  (B.2).  For  determination  of  the  load  and 
bending  moment  required  to  Join  the  gussets  and  baffle,  the 
disparity  AR^  from  Eq.  (B.2)  is  equated  to  the  sum  of  AR^ * , 
AR2  and  AR3  as  given  by  Eqs.  (B.4),  (B.6)  and  (B.9)  respec- 

t i vel y ,  or 
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(B.17) 


ARX  =  AR'  +  AR2  +  AR3 

This  allows  evaluation  of  W,  M.  M' ,  M  and  VV  so 

7  7  7  o  o 

that  the  combined  stresses  in  the  brake  can  be  determined. 

a )  For  the  Conical  Baffle: 

With 

P  =  5000  psi 

R^  =  4.16  in.  (mean  radius  of  1"  section 

at  outer  diameter  of  baffle) 

t  ^  =  0.98  in. 

E  =  30  x  106  psi 

Y  =  62° 

from  Eqs.  (B.2)  and  (B.4) 

ARX  =  0.0062  in.;  ARX  *  =  6.45  x  10_1°  W 

b )  For  the  Gussets 

With 

A  =  2.5  in. 

I.  =  0.5  in.4 

2 

A  =  1.5  in. 

Z  a 

G  =  llx  10  psi 
from  Eq.  (B.6) 

AR2  =  5.28  x  10'  7  W 

c )  For  the  Mounting  Rine 

With 

9-  -  2.5  in. 

R3  =  3.04  in. 

t3  =  0.70  in. 

V  -  0.25 

X  =  0.891 

D  =  9.10  x  105 
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AR3  =  25.0  x  10'  7  W 

Substituting  these  values  into  Eq.  (B.17)  results  in 

W  =  2260  lb 

ARj  *  =  14.7  x  10'7  in. 

AR2  =  0.0012  in. 

AR^  =  0.0050  in. 

This  indicates  radial  deflection  of  the  baffles  is  not 
appreciably  restricted  by  the  gussets  and  mounting  ring. 

B.4.2  Evaluation  of  Combined  Stresses 

a  )  Conical  Baffle 

From  Eqs.  (B.l)  (B.3)  and  (B.5) 

51  =  45,000  psi 

S  •  =  8,000  psi 

Maximum  circumferential  tensile  stress  at  inner  surface  of 
baffle,  Sx  +  S  •  =  53,000  psi. 

b )  Gussets 

From  Eqs.  (B.7)  and  (B.8) 

S  £  =  45,700  psi 

Sb2  =  11,400  psi 

Maximum  tensile  stress  at  Junction  between  gusset  and  mounting 
ring  inner  surface  becomes 

52  *  Sb2  =  57>100  Psi 
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Tensile  stress  in  same  section,  upper  surface,  becomes 

S2  '  Sb2  =  34>300  P®1 

c )  Mounting  Ring 

From  Eqs.  (B.15)  and  (B.10), 

Sfa3  =  10,500  psi 

S^3  =  19,400  psi 

These  maximum  stresses  occur  at  different  locations  in  ring 
and  cannot  be  combined. 

d )  Tube  Extension 

Using  Eqs.  (B.l)  and  (B.2)  with  subscripts  4, 

R.  =  2.39  in. 

4 

t  .  =  0.33  in. 

4 

Y  =  0° 

S4  =  30,200  psi 

AR„  =  0.0029  in. 

4 


a  nd 
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PRESSURE  AND  FLOW  FIELD 

EXTERIOR  TO  A  GUN  BLAST 

By 

A.  C.  Eringen  and  J.  D.  Ingram 

C .  1  PHYSICAL  BASIS  OF  THE  PROBLEM 

The  health  and  safety  of  a  gun  crew  requires  that  the 
acoustic  noise  level  produced  by  a  gun  blast  not  exceed  a  pre¬ 
determined  value.  To  a  first  approximation,  the  noise  level 
is  proportional  to  the  pressure  rise.  The  Human  Engineering 
Laboratories (l  )  recommend  that  no  persons  should  be  subjected 
to  an  overpressure  exceeding  7  psi.  Thus  the  problem  of 
determination  of  the  pressure  fields  and  flow  exterior  to  the 
gun,  following  the  blast,  must  be  solved.  The  basic  problem 
may  be  divided  into  two  parts:  (a)  the  flow  and  pressure  field 

interior  to  the  gun  tube;  and  (b)  the  flow  and  pressure  field 
exterior  to  the  gun. 

The  solution  of  the  problem  (a)  supplies  us  with  the 
initial  conditions  necessary  for  the  problem  (b)  as  may  be 
seen  from  the  following  consideration:  Upon  ignition,  the 
shell  moves  as  a  piston  in  the  gun  tube,  thus  pushing  the  air. 
This  air  is  compressed  so  that  a  shock  wave  proceeds  toward 
the  muzzle.  The  flow  condition  between  the  shock  wave  and  the 
shell  and  the  waves  generated  and  reflected  between  the  two 
must  be  studied  in  order  to  know  the  initial  condition  at  the 
outlet  of  the  muzzle.  This  problem  has  been  studied  extensive¬ 
ly,  (2)  and  (3),  in  the  form  of  various  simplified  mathematical 
models.  The  full-fledged  problem  is  highly  complicated,  but 
not  impossible. 
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When  the  shock  first  reaches  the  muzzle,  suppose  that  the 
flow  and  pressure  fields  at  the  muzzle  are  known: 


Figure  Cl 

That  is,  suppose  that  the  pressure  pQ,  density  pQ,  and 

velocity  IT  at  the  origin  of  our  rectangular  coordinates 
o 

x,  y,  z  are  determined  through  problem  (a).  We  select  the 
origin  of  time  when  shock  waves  reach  the  muzzle  of  the  gun* 

Given  the  initial  state  p^,  pQ,  u^  the  problem  (b)  is 
to  determine  p,  p,  tT  at  any  space  point  as  a  function  of  time* 

It  is, theref ore ,  clear  that  the  problem  (b)  cannot  be 
solved  independently  of  problem  (a),  and  that  they  are  coupled 
through  the  muzzle  shock*  In  the  treatment  of  problem  (b)  we 
assume  that  p^,  pQ  and  uq  at  t  =  0  are  given*  The  exterior 
flow  field  now  requires  the  solution  of  the  basic  equations  of 
gas  dynamics  subject  to  some  initial  and  boundary  conditions. 


- <QE  A> - 
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The  Initial  Conditions  require  that  we  prescribe  p,  p,  IT 
at  time  t  =  0  throughout  the  region  exterior  to  the  gun  tube. 

The  Boundary  Conditions  require  that  on  the  lateral  surface 
of  the  gun  and  on  any  other  rigid  surface  (such  as  the  ground) 
the  normal  component  of  the  velocity  vanishes  and  at  infinity 
the  air  be  at  rest  at  constant  pressure  and  density. 

The  blast  waves  from  a  spherical  explosion  have  been 

treated  by  Taylor  (4)  and  others  (5),  (2).  However,  in  the 

treatment  of  this  problem,  the  so-called  self  similar  motion 

forms  have  been  employed.  It  is  known  that  such  solutions 

2 

break  down  at  pressure  levels  of  4kg/in.  .  Moreover,  the 
directional  effect  of  the  gun  blast  is  not  included  in  station¬ 
ary  point  blast  waves.  Clearly  the  pressure  at  points  which 
are  mirror  images  of  each  other  with  respect  to  z  —  0  plane 
are  different.  Therefore,  spherical  blast  waves  are  not 
acceptable  from  the  point  of  view  of  their  lack  of  directionality. 

The  acoustic  approximation  treated  in  (6)  on  the  other 
hand  considers  only  motions  due  to  density  and  heat  addition 
inside  a  spatial  region.  The  present  problem,  however,  is 
basically  an  i ni t ia 1 - va 1 ue  problem,  requiring  a  different 
treatment.  In  order  to  cast  light  on  the  directionality  of 
the  problem  and  provide  some  intuitive  basis  for  the  most  exact 
problem  of  finite  amplitude  waves,  we  give  here  a  treatment  of 
the  limiting  case  of  acoustical  approximation. 

C  .  2  ACOUSTICAL  PROBLEM 

The  linearized  equations  of  an  axisyrametric  acoustic  field 
expressed  in  cylindrical  coordinates  (r,  z)  are  (see  Fig.  Cl): 


- <Q  E  A> - 
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0 


l£ 

dt 


+ 


a 


M  .  i_  il£ 
at  Pa  3r 


0 


.  ,L_  &JB. 

at  p  dz 

a 

2 

P  =  ca  P 


0 


(C2.1  ) 


Here  the  first  equation  is  the  equation  of  conservation  of  mass; 

the  next  two  are  those  of  the  balance  of  momentum  and  the  last 

equation  is  the  equation  of  state.  The  quantities  p,  p  and 

(u,  w)  are  respectively  the  pressure,  the  density  and  the 

velocity  components.  The  ambient  density  and  the  constant 

velocity  of  sound  are  denoted  by  p  and  c  respectively. 

a  a 

In  vector  notation  these  equations  may  be  written  as 


+  p  divtT  = 
dt  r  a 


3u  la_ 

dt  +  ~  grad  P 

a 


=  0 


(C2.2) 


P  =  ca  P 


Initial  Conditions  require  that  at  time  t  =  0  we  prescribe 
the  velocity  and  pressure  fields,  i.e.. 


u (r ,  z  , 

0) 

- 

u  (r ,  z  ) 
o  7 

w  (  r  ,  z  , 

0) 

= 

w  (  r  ,  z  ) 
o  7 

CC2. 3  ) 

p(r,  z. 

0) 

= 

PQ ( r ,  z)  =  ca2  pCr,  z. 

0) 

2  , 

=  c  p  ( r , 

a  r  o  7 

w  a  nd 
o 

Po 

are 

1  given  functions. 

_ yrrprs 

N  C  E 

E  N 

6  1 

yu  l  r\  p 
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The  effect  of  the  boundaries  of  gun  tube  will  be  neglected. 
The  problem  is,  therefore,  given  an  initial  velocity  and  pressure 
(or  density)  distribution,  to  determine  the  resulting  acoustical 
field.  The  solution  is,  further,  subject  to  vanishing  velocity 
field,  constant  atmospheric  pressure,  and  density,  at  infinity. 

C .  3  THE  SOLUTION  OF  THE  ACOUSTIC  PROBLEM 

For  the  treatment  of  the  i ni t ia 1 -va 1 ue  problem  at  hand, 
the  most  appropriate  tool  is  the  Laplace  transform  as  given 
in  (7) 

F 

F(r,  z,  s)  =  /  e  8t  F(r,  z,  t)  dt  (C3.1) 

o 

which  has  the  inversion  integral 

rT  +  iP 

F(r,  z ,  t)  =  I  ®8t  F(r,  z,  s)  ds  (C3.2) 

r  -  i? 

where  y  a  nd  P  are  real.  For  the  region  and  conditions  of 
validity  of  (C3.1)  and  (C3.2),  see  (7). 

Applying  the  Laplace  transform  to  (C2.2)  we  get 

sp+p^divu  =  p^ 

2 

c 

su  +  ~ —  grad  p  =  XT  (C3.3) 

pa 

2  - 

P  =  ca  P 

where 
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oo 


p  =  p  (  r  ,  z  ,  s  )  = 


u  = 


=  u  ( r  ,  z  ,  s)  = 


-at 

p  (  r  ,  z  ,  t)e  dt 


-*■  v  -st 

u  (  r  ,  z  ,  t)e  dt 


(C3.4 ) 


etc.,  and 

tT  =  u  e  +  w  "cT 


"ft  =  u  7T  +  w  e* 
o  or  o  z 


(C3. 5  ) 


are  the  velocity  vectors  and  e^  and  are  the  unit  vectors 

in  r  and  z  -  directions. 

From  the  last  two  of  (C3.3)  we  solve  for  TT  : 


tT  =  ( tT  -  -  grad  p) 

s  o  p  K 


(C3.6) 


Substituting  this  into  the  first  of  (C3.3)  we  get 


V2  p  -  (s/c  )2  p  =  F(r,  z,  s) 
a 


where 


v2  S  =  iiii  +  1  fi  + 

-  2  r  dr  $z2 


f 

dr 


(C3. 7) 


(C 3 . 8  ) 


F(r,  z ,  s )  =  P  div  IT  (r,  z)  -  (s/c  )  p  (r,  z) 

7  7  a  o  7  a  o 


The  problem  is  now  reduced  to  solving  the  inhomogeneous 
Helmholtz  equation  (C3.7)  for  the  infinite  domain.  This  may 
be  achieved  by  determining  the  Green’s  function 
G(r,  6 ,  z  {  JE  |  (p ,  £  ,  t)  whose  Laplace  transform  satisfies 
(C3.7)  with  F  replaced  by  the  three  dimensional  delta  func¬ 
tion,  i  .  e  .  , 


V2  G(r,  0,  z;  £  ,  <P  9  £  ,  s)  “  (s2/c^2)  G  =  S(r-|,  z-£,  B-<p) 


2  y  2,  - 
a 
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The  solution  to  this  equation  in  three  dimensions  is  well-known 
It  is  given  by 


-kR 

G  =  — - 

47TR 

R2  =  r2  -  2r£  cos(0-^)  +  £2  +  (z-£)2 

k  =  (s/c  ) 
a 

The  solution  to  equation  (C3.7)  may  now  be  expressed  as 


P  = 


s)  G(r,  0,  z;  g  9  f  £  f  s)£d|d^?d£* 


Since  F  does  not  depend  on  <p  we  may  integrate  on  <p 
To  this  end  we  make  use  of  the  Sommerfeld  representation 
point  source  at  (£,  i.e., 

jx 


-kR 


JQ(rp)  e 


-  IWIV  Z 


dr 


where 


■')2  =  t2  +  *2 

2  2  2 
p  =  r  -  2r|  cos(0-^O  +  | 


The  addition  formula  for  Bessel  functions  is 


jo<tP) 


00 


n=0 


£  n  J  C^f)  Jn(2T)  008  n(0^^} 


1 

2 


n  =  0 

n  =  I  p  2  f 


Thus 


- <Q  E  A> - 
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(C3.9) 

(C3.10) 

directly, 
of  the 

(C3.ll ) 


(C3.12) 


I  N  C. 


-kR 


oO 


oo 


€n  Jn(r#)  Jn(rr)  cos  n(0-p> 


-y\z-ci 


n=0 


Integrating  on  <p  from  0  to  2tt,  only  the  term  for  n=0  con¬ 
tributes  so  that  we  obtain 


G* 


crp  Jo(rr> 


The  function  F(r,  z,  s)  given  by  (C3.8)  is  composed  of 
two  terras  one  of  which  is  linear  in  s  and  the  other  indepen¬ 
dent  of  s.  Consider  now  the  following  Laplace  transform: 


dt 


-P(0)  +  sP 


(C3  # 12a ) 


In  our  case,  the  initial  pressure  and  velocity  are  prescribed 
inside  a  closed  region,  outside  of  which  they  both  are  zero* 

Hence,  we  may  invert  the  Green's  function  as  if  it  were  inde¬ 
pendent  of  s  and  for  the  terifi  linear  in  s  simply  take  the 
time  derivative  of  the  resulting  expression*  Thus  if  G*  is 
the  Laplace  transform  of  the  Green's  function  for  the  axisymmetric 
inhomogeneous  Helmholtz  equation  which  is  independent  of  s,  the 

response  to  an  initial  pressure  p  and  velocity  tT  is 

o  o 

P<r,  *.  t)  =  JJ  <paV  •  G*>{dfdf  -  fj 

ca 


where 

G*(r,  z;  |  t)  =  L_1  jc*<r,  z,|  s)f 

is  the  Inverse  Laplace  transform  of  5*. 


(C3.13  ) 
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ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


C  -  8 


Returning  now  to  the  axially  symmetric  Green's  function,  we 
found  that 


oo 


6*  =  V  Jo(Z'()  Jo(rr) 


-V  |  z-C\  rdg 
’  V 


Now,  apply  the  inverse  Laplace  transform  to  this  equation 

°*  =  J  3o(tl)  Jo(2'r)  *  |a'<rl  *?*} 


*  <‘5>  V^)  V”>  L_1 


-*V|  «-<C| 


■)  z ar 


(C3.14 ) 


^00 


*  ’<?> 


J  rcP  Vz'r)  vW*?'2  ■  >  roe*  <=< 


=  0  >  0<cat<lz~^  I 


We  may  evaluate  this  Integral  exactly  (8) 


G*  = 


4ttA 


(C3.15 ) 


»  r»  ]f 


if 

otherwise . 


2  *  2  ,  ^X2 

c  t  (z-s) 

a 


Here 


* 

G-0 


make  up  a  triangle,  and 
A  is  defined  as  the  area  of  the  triangle  whose 


/ 


2  2  2 

sides  are  £  ,  r,  and  ^c  t  “  (z-^)  •  As  for  the  conditions 

for  these  to  compose  a  triangle,  first,  they  must  all  be  positive 
second,  they  must  satisfy  the  following  inequalities: 


+  r  + 


i 


2.2 
c  t 
a 


-  U-gY 


i- 


r  + 


ca2t2 


-  (z-£T 


(a) 


(b  ) 


(C3.16 ) 
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(c) 


<r+f)  -  /ca2t2  -  (z-£)2  »  0 


(d) 


These  then  are  the  defining  equations  of  the  region  of  definition 
of  the  Green's  function.  The  region  R  defined  by  (C3.16  b, 
c,  d)  is  the  inside  of  a  circle  in  the  |  plane  having  its 
center  at  (r,  z)  and  with  radius  c  t.  If  the  non-zero 
region  of  initial  pressure  and  velocity  is  denoted  by  S,  then. 


for  the  response  to  an  arbitrary  p 


o 


r  n  s 


Po(if^)  dt^  I  d!d£  (C  3 . 1  7  ) 


i 


2 


a  Rf)S 


where  RAS  signifies  the  region  common  to  both  R  and  S.  (See 
Fig.  C2  ) 


Figure  C 2 
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The  explicit  form  of  (C3.15)  is 


-  (c  /t r) 

: = . =  - r . . .  — 

~  "  <r-f)2  \J  Cr+|)2  +  (z-£)2  -  ca2t2 

(C3. 18 ) 

Equa t i on  (C3 *  1 7 )  gives  the  pressure  field  for  an  arbitrary  initial 
pressure  and  velocity*  The  velocity  field  may  now  be  determined 
from  (C3.6).  This  is  omitted  here* 

C  *  4  Special  Cases 

C  *  4  „  1  Case  I  -  Specification  of  Initial  Velocity  Distribution 

Consider  an  axial  initial  velocity  distribution  over  the 
cylindrical  region  -  oo  <  z  <  zq,  0  <  r  <  a,  i.e., 

u“  =  t  fv  (l-U(z))  +  V  (1-z/z  )  (U(z )  -  U  (z-z  )]  ,  0<r<a  (C4.1) 

o  z  L  o  o  o  o  J  ' 

=  0  otherwise 

where  U(x)  is  the  unit  step  function  defined  by 

U  ( x )  =  0  x<0 

U  ( x )  =  1  x>0 

The  velocity  distribution  (C4.1)  is  shown  in  Figure  C3. 


Figure  C3 
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From  (C4.1)  it  follows  that 


— ►  v  _ 

V  °  If*  =  -  — -  u(z)  -  u (z-z  )] 

O  2  l  O  J 


The  region  S  then  becomes 


0  < 


0  <  I  ^ 


From  (C3017)  and  (C4.1)  it  follows  that 


P  = 


P  v 

ra  o 


j 


ROS 


G*fd4d£ 


(C4 . 2  ) 


(C  4 . 3  ) 


This  integral  may  be  integrated  with  respect  to  £  immediately, 


where 


P  = 


c  p  V 
a  r  a  o 

TTZ 


ta  n 


-1 


<r 


t2  -  a2 

2  -2 

B  -  l 


d£ 


(C4.4) 


max  t  (£ ) 
ROS 


A  = 


B  = 


yTV  -  (z-<f)2  -  r 


/ 


caV  -  (z-<)2  +  r 


(C  4  #  5  ) 


and  £  and  C  are  determined  by  maximum  and  minimum  £  possible 

in  ROS,  (Figure  C4)0  The  term  corresponding  to  the  lower  limit 

2  2 

in  the  £  integration  vanished  since  J.  =  A  •  For  the  case  of 


( 


c_ 2t2  -  (z-£)^  +r  <  a  we  have  ^=B#  Thus 
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Figure  C 4 


h  oweve  r , 

B>a  then 

we  have  £ 

:  =  a 

a  nd 

c 

p  c  V 
ra  a  o 

1  tan  ^ 

J£ 

U2 

Cl 

< 

P  = 

WZo 

u2 

2 

-  a 

C4. 7  ) 


This  i s  the  case  for  r>a.  In  order  to  evaluate  this  integral, 
we  integrate  once  by  parts  to  obtain 


P 


p  c  V 

ra  a  o  .  -1 

(z-£)  tan 


7TZ 


p  C  V 
ra  a  o 

7TZ 
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Now 


B  -1  / a2-A2 
*  V  B2-a2 


-<z-<p  [  ca2t2-(Z-<)2-  r2+a2] 


2[ca2t2-(z-C)2]  ^ca2t2-(z-^)2-(r-a)2  ^(r+a  >2+ (z-^-c^t2 


(C4.8) 


and  therefore. 


facaVo  .  -1  /a2-A2 

P  '  -  ^  (Z^>  1/  15 2 

o  V  B  -a 


c=  i 


<r=  c 


p  c  V  I  <Z-^)2  [  c  2t2-(z-^)2  -  r2+a2l  d< 

a  a  o  | _  La  ■J  ^ 


27TZ 


ca2t2-(z-^)2]  yCa2t2-(z-^)2-(r-a)2^(r+a)2+(z-^)2ca2t2 


(C4.9) 


Let  us  examine  the  Integral  in  detail.  Let 
z-C 

_  I 


ft2  [ctV-^rV]dA 


U  [f-%2'2]  /c.V-f-Cr-a)2  W*<  r->2-c.2t: 


(C4  •  1  0  ) 


This  integral  may  be  expressed  in  terms  of  elliptic  functions 
We  will  make  use  of  the  following  expressions  (9)* 

[  a2  dA 


-  C2)  (b2  -  ^2) 


=  gt>  E(^, 


k) 


(C4. 11 ) 
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/ 


q2  -Pj)dA 


b2-P 


(^2-P)  / (b2-^2)(^2-c2) 


—  *  ^  [(a2-a12)7r(^,a2,k)+a12F(^,k)] 


(b  -P)  a 


where 


1  .  2  _  2  2 ,  .2 

g  =  —  k  =  (b  -c  )/b 


2  ,.2  a 2  w  /u2  2  N 

sn  u  =  (b  - /.  )/(b  -c  ) 


(C4.12) 


=  sin 


-1  /b2-v2 

b2-c2 


2  _  (b2-c2) 


2  _  (b2-c2) 

12  ®  “  2 

(b  -P^  (b  -P) 


F  (^,  k),  E  (j ir,  k)  and  7 r(^f,  a  ,  k)  are  the  elliptic  integrals 
of  the  first,  second  and  third  kind  respectively,  and  are  defined 
by 

f 


<.f,  k>  = 


d© 


v*- 


2  2 
k  sin  © 


E(f,  k) 


-  /  f~. 


2  2 

k  sin  6  d0 


(C4.13  ) 


a2,  k)  =  J 


dJL 


(1  -  a2sin20 


)  iA- 


k2sin20 


Applying  these  to  equation  (C4.9)  we  obtain 


P  = 


p  c  v  p  - 
a  a  —  ^(z-4  )  tan 


TTZ 


/a2-A2 

§2-a2 


-  (z 
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p  c  v 
_  a  a 

2t rz 


2  / 


2.2  ,  \ 2 

ca  t  -  (  r-a ) 


,  k)  -  E(^,  k) 


2  2 

2  2.  c  t 


+  p7fc^ilF<£  °2’ k> "  a2>  k)l 

+  L,*]{r<,t’ k)  -T<h  k>i] 


(C  4  •  1  4  ) 


where  the  quantities  having  bars  above  or  below  denote  the  corres¬ 
ponding  quantities  where  £  is  replaced  by  <£"  and  £.  respectively, 


C • 4 « 2  Case  II  -  Specification  of  Initial  Pressure  Distribution 

As  shown  in  equation  (C3.12a)  we  need  only  multiply 
2 

(C4.14)  by  P  /c  p  V  )  and  differentiate  with  respect  to 

o  a  1  a  o 

time  to  obtain  the  response  to  an  initial  pressure  distribution 

Po  =  "  U(z-zq)J  r  <  a 

=  0  r  >  a 


(C4. 15 ) 


Thus  we  obtain 


P  = 


-P 


7TC 


i-f 

st  L 


(z-Otan 


-1  /a2-A2 
VS2-a2 


-  (z-^)tan  1  ,/a2~A2 
B  -a 


i 

J 


k) 


+  _ (r-aj. 


[..V-c-o2] 


>F(^,  k)  -  F  (if[,  k)i 


(C4. 16  ) 
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It  l 8  possible  to  write  the  indicated  partial  differentiation  in 
explicit  forms*  However,  the  expressions  so  obtained  are  cumber¬ 
some  and  for  the  present  numerical  calculations  were  not  needed. 

C . 5  ASYMPTOTIC  SOLUTIONS 

The  foregoing  exact  solution  of  the  acoustic  problem  for  a 
gun  blast  is  expressed  in  terms  of  complicated  elliptic  functions 
which  are  not  readily  amenable  Xo  calculations.  A  limited  amount 
of  calculations,  based  on  the  exact  solution  (C4.14),  are  carried 
out  on  a  digital  computer  and  plotted  in  Figures  C6  to  C8. 
However,  in  order  to  gain  a  greater  Insight  into  the  nature  of 
the  problem,  it  is  felt  that  an  asymptotic  approach  may  be  use¬ 
ful.  The  far  field  approximation  is  the  basic  theme  of  this 
article.  To  this  end  we  explore  the  possibility  of  finding  a 
Fourier  integral  representation  to  the  solution  of  our  basic 
equation  (C3.7).  We  divide  the  space  into  two  regions  r£a  and 
r>a ,  and  find  the  solutions  of  (C3.7)  subject  to 

F  =  0  for  r>a  (C5.1) 


and  that  the  pressure  p  and  its  partial  derivative  3p/gr  are 
continuous  at  r=a. 


The  Fourier  transform  <]5(s)  of  a  function  $(z)  and  its 


inverse  are  respectively  given  by 

oo 


$(£)  =  ~  I  $(2)ei^Zd2 


-  CD 
00 


(C5.2  ) 


$(z)  =  1  <£(Oe~lCzd< 


-oo 
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We  now  try  to  find  a  Fourier  integral  representation  for  the  solu¬ 
tion  of  (C3.7)  in  the  forms 

OO 


P  =  |<<)  K  (Vr)  e  ls*d£  r^a 


(C5.3) 


-oo 

oo 


00 


p  = 


yS)  Io(^r)  e'^Zd^ 


f  (s  )e~*^Zd£ 


2  2 

k  +£ 


-oo 


—00 


r<a 

(C5.4) 


where  £(*f)  and  1J  (£ )  are  two  functions  to  be  determined, 


I  (Vr)  and  Ko(Vr)  are  the  modified  Bessel  functions  of  zeroth 


order  and 


V- 


k  =  s/c 


a 


oo 


(C5.5) 


Us*C)  =  ^  /  F(s,z)e‘^Zd£ 


-oo 


Here  we  have  assumed  that  F  is  independent  of  r  for  r<a. 


The  continuity  of  the  pressure  p  and  its  normal  gradient 
3p/^r  at  r=a  and  t>0  requires  that 


i  <<)  K0<Va)  =  rj(£)  Io(Va)  -  f(s,^)// 

|(4D  Kx(Va)  =  -77(C)  ^<7*) 

Solving  these  two  equations  for  £  and  ^  we  find 


f(C)[Ko(Va)  Ij^CVa)  +  J^CVa)  I  q  CV  a  )J  =  -  I^a) 


or 


|(C)  =  -  ^  f(s,C)  1 !  (*Va  ) 


(C5 . 6 ) 
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where  we  have  employed  the  Wronskian  relation 


K  (x)  I.  (x)  +  K  (x)  I  (x)  =  - 

O  1  1  O  X 

Sirailarl y 

rj(£)  =  f(s,^)  Kx(Va)  (C5.7) 

The  final  solution  is,  therefore,  found  to  be 

CD 

P  =  -L'1  |  J  a  f(Sv^J>  Ix(ya)  Kq  CVr  )e~^^"Zd^" 

-  CD 
00 

„  .  l-1  | J  [a  .oCVD  K^a)  - 

-  oo 

It  is  these  integrals  which  we  will  use  for  the  asymptotic  devel¬ 
opment.  The  pressure  field  will  be  determined  at  points  far  from 
the  source,  i.e,  for 


r>a 


(C5.8) 


]  .-«■*] 


r<a 


R 


✓ 


(r-a  )4 


2 

z 


»  1 


In  general  £  and  s  are  complex  numbers  such  that  the  integrals 
appearing  in  (C5.8)  are  exponentially  dominated.  Asymptotic 
development  of  integrals  of  the  form 
oo 


I 


K  Hr) 
o _ 

V 


i1<.y&)e'i€z  d£ 


(C5.9) 


-  CD 


is  obtained  for  the  region  r^a .  We  rewrite  this  integral  in  a 
form  convenient  for  our  purpose. 

oo 

[eN(r-a)  Ko(Vr>  IX<V«)J  *£•  (C5.10) 

-GO 
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The  quantity  inside  of  the  bracket  in  the  integral  is  a  slowly 
varying,  well  behaved  function®  Replacing  the  Bessel  functions 
by  their  asymptotic  expansions  we  obtain 

CD 


I 


-V(r-a) 

e 


-  OO 


- - -  /  e  -  V<r-a  )  -  i<-z  d£  r>& 

2  \ZrT  J  y 


(C5* 11  ) 


-  oo 

The  method  of  steepest  descent  may  now  be  applied  to  this  integral 
to  obtain  its  form  for  large  R.  We  will  find,  however,  that  the 
saddle  point  lies  close  to  the  pole  introduced  by  the  use  of  the 
asymptotic  expressions  for  the  Bessel  functions*  This  pole  is  not 
a  part  of  the  solution  and  may,  therefore,  be  removed  by  the  use 
of  a  Laurent  series  expansion*  First  we  examine  the  position  of 
the  saddle  point*  The  exponential  portion  of  the  above  integrand 
i  s 


E  =  -  V  (r-a  )  -  i£z 

At  the  saddle  point  we  have 


af  =  f  <r-a>  -1*  =  ° 

The  solution  of  this  equation  gives  the  saddle  point  g 


L  - 


-ikz 


\J77- 


>.2  2 
a  )  +z 


- 1  kz 

R 


(C5. 12  ) 


At  the  saddle  point,  takes  the  value 

*0  =  XiT-L.a-1  (C5.13) 

ro  R 

2 

Since  the  pole  corresponding  to  1/y  is  a*  =  “ik»  the  saddle 
point  lies  between  the  origin  and  the  pole  as  shown  In  Fig.  C5  . 
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Fig.  C 5  -  COMPLEX  -  PLANE 


As  stated  above,  we  now  remove  the  singularity  at  the  pole  by 

2  2 

expanding  l/(£  +  k  )  in  a  Laurent  series,  i.e., 


_ 1 _  _  _ izl/21 _  .  _i_  .  LL±  ik) 

2  2“  2  ^  •  •  * 
£  +  kz  +  ik)  ik  4  k  8ikj 


Ignoring  the  terms  in  when  n>0,  we  obtain 

oo 

JL  =  1  =  (-1/2)  +  1  (i}n 

V2  <2+k2  (£*  +  ik)  ik  4k2  ^  2 

n=0 


2ik( £  +  ik)  2k2 

Integral  (C5.ll),  with  this  substitution  and  dropping  the  residue 
contribution  due  to  pole,  gives 


I 


r 

r  k3  <*-<o  >1 

t.  "kR 

‘Ik(r-a)2  2  J 
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where  we  have  expanded  the  exponent  in  Taylor  series  about  the 
saddle  point  and  retained  only  terms  up  to  second  order#  Making 
the  substitution 


Co  =  ft  x 

d£  =  Jk  dx 

and  extending  our  path  of  integration  to  infinity  in  both  direc¬ 
tions  along  the  path  of  steepest  descent  we  find 


I 


4 


-kR 

fa  k3/2 


-R 


(r-a) 


_  x 

2  2 


dx 


-oo 


^  Air  *  l 

4K3''2  „3/2 


-kR 

e 


(C5.14) 


The  asymptotic  representation  of  the  pressure  field  is,  therefore^ 
given  by 


p (r , z , t  )  = 


<r-a  > 

4  /£  7T  C 


OO 


-kR 


o-  kl/2  r3/2 


d£ 


-oo 


where 


R  = 


z)2  +  (r  -  a  )2 


-1 


Upon  exachanging  the  operator  L  and  the  integral  we  find 

R 


oo 


-S“ 


p(r,z,t)  = 


JUi=!L 


/ 


4  ,/277-C 


P0(O 

3/2 
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-oo 


1/2 


d<T 


ORDNANCE  ENGINEERING  ASSOCIATES,  INC. 


C  -2  2 


«■  n T  r  *  ’  fTF”vr  *•’  nr  *■  *wwiu  \  f 
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U ( t-R/c  ) 
_ _ a  __ 

^7rTt-R/ca  ) 


where  U(t-R/c  )  is  the  unit  step  function.  Consequently  the 
a 

asymptotic  approximation  of  pressure  field  takes  the  final  form, 

oo 


P 


(r-a) 

/2  4tt 


-  00 


p  (/“)  U  ( t-R/c  ) 
_  _  _  _  _  a  _ 

R3/2  \/t-R/ca 


(C5. 15  ) 


We  may  now  consider  some  special  cases  of  pressure 
fields.  Below  is  presented  the  case  of  uniform  initial  pres¬ 
sure  over  a  finite  region  CO,  z  ) ,  i.e., 

Po(£)  =  PQ[u(^)  -  U(<^  -  2q)J  (C5.10) 


where  P  is  a  constant  and  U(x)  is  the  unit  step  function. 

o 

Substitution  of  (C5.16)  into  (C5.15)  gives 


oo 


-oo 


The  first  pressure  arrives  when 


2  2  ,  .2  ^  ,  n2 

c  t  =  (r-a  )  +  (z-z  ) 

a  o 

Theref  ore 
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from  then  on 


P  = 


s/2  4-tr  '  *  J  \f 


'cQt-R  R 


3/2  > 


c  t  >. 

A  f 


\[77- 


n2  2 

a  )  +z 


We  can  see  quickly  that  the  maximum  of  this  integral  occurs  at 
2  2  2  2 

c  t  =  (r-a)  +  z  and  is  given  by 

a 


p  (  r  ,  z  )  = 

v  7 


max 


\[2  4tt 


$ 


(r-a  ) 


—  3=^ 

o  J(  r  -a  )2  +  (z-£*)2  JJ(r-a  )2  +  z2  -  ^(r- 


-a)2+(z-£)2 


Obviously  most  of  the  contribution  to  this  integral  occurs  at 
£  =  0 ,  he  nee 


p  (  r  ,  z  ) 

max  7 


(  r -a  ) 


» 2  2 
a  )  +  z 


1 


This  expression  is  valid  for  z^z^*  Similar  ones  may  be  obtained 
for  0<z<zq  and  z<0.  These  asymptotic  solutions  do  not  have 
a  well  defined  domain  of  validity*  However,  we  may  make  the 
qualitative  statement  that  they  break  down  for  z>>r*  This  is 
caused  by  the  use  of  the  asymptotic  form  of  the  Bessel  functions* 
An  approximate  region  of  validity  is  indicated  by  the  dotted 
lines  in  Fig.  Cll.  Plots  of  the  isobars  for  the  peak  in  the 
pressure  distribution  are  given  in  Fig.  Cll* 


C . 6  CONCLUS IONS  AND  SUGGESTIONS 

The  results  of  the  numerical  work  done  on  this  problem  are 
given  in  Figs.  C6  to  C12*  Figures  C6  and  C7  show  the  waveforms 
obtained  from  equations  (C4.14)  as  a  function  of  the  radial  and 
axial  coordinates.  Note  that  the  decay  in  the  z  direction  is 
more  rapid  than  in  the  radial  direction.  This  is  a  manifestation 
of  the  development  of  cylindrical  waves.  This  phenomenon  is 
shown  clearly  in  Fig.  C8  in  which  the  isobars  for  the  peak  in  the 
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waveform  are  given.  Figures  C9  and  CIO  are  the  waveforms  ob¬ 
tained  from  the  asymptotic  development*  It  should  be  noted 
that  these  are  infinite  in  extent  although  they  decay  in  time. 

The  isobars  for  the  asymptotic  development  are  shown  in  Fig. 

Clio  Here  the  region  of  doubt  for  these  are  indicated  by  dashed 
lines.  The  decay  along  the  z  axis  for  specified  initial 
velocity  distribution  is  given  in  Fig.  C12. 

The  analytical  work  presented  here  can  be  applied  to  other 
types  of  initial  velocity  and  pressure  distributions.  Further 
numerical  work  will  no  doubt  cast  better  light  into  the  nature 
of  the  problem.  The  acoustical  problem  of  gun  blast  with 
muzzle  brake  can  also  be  formulated  and  solved  with  the  technique 
employed  here.  The  effects  of  the  presence  of  the  gun  and  the 
reflection  of  the  wave  from  the  ground  need  novel  considerations 
not  included  in  the  present  analysis. 

Finally  a  work  or  two  may  not  be  out  of  order  in  regard  to 
the  nonlinear  problem  involving  strong  shocks.  From  the  above 
solution  of  the  acoustical  problem  it  is  clear  that  the  exact 
solution  of  the  correct  nonlinear  problem,  presently,  is  out  of 
our  reach.  However  an  approximation  in  the  form  of  spherical 
shock  waves  emanating  from  a  point  explosion  with  a  given  energy 
can  be  treated.  While  certain  solutions  of  this  problem  are 
found  in  the  literature,  (2),  (3),  (4),  they  are  not  directly 

applicable  to  the  problem  at  hand  since  no  directionality  is 
inherent  in  these  solutions.  The  final  solution  of  the  non¬ 
linear  problem  requires  a  skillful  combination  of  analytical  and 
computer  work. 
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RELATIVE  PRESSURE 


c  t  -  |/(r-a)2  +  (Z-Z  )Z 

t*  -  a  ■  r  ,  ,  — - -9 — 

/(r+a)2  +  (Z)2' -  /(r-a)2  +  (Z-Zq)2 


Fig.  C 6  -  EXACT  WAVEFORMS 

p  _  (Pleasure  at  r  and  z) 
r  (Pressure  at  r=5',  z=5’) 
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RELATIVE  PRESSURE 


jr  -1 


C  t  -  \/(r-a)2  +  (Z-Z  )2 

t*  =  -  V -  o 

/(r+a)2  +  (Z)2  -  |/(r-a  )2  +  (Z-Zq)2 


Fig.  C  7  -  EXACT  WAVEFORMS 

p  _  (Pressure  at  r  and  z) 
r  (Pressure  at  r=5 ' ,  z-5 *  ) 
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Fig.  C  8  -  PEA  K  _PRE.S_S.URE.,.  jS  Op  A  R  S 
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Fig.  C  9  -  ASYMPTOTIC  WAVEFORMS 

p  _  (Pressure  at  r  and  z) 

r  "  (Pressure  at  r=5 * ,  z'=5') 
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Fig.  CIO  -  ASYMPTOTIC  WAVEFORMS 

p  _  (Pressure  at  r  and  z ) 
r  (Pressure  at  r=5 * ,  4=5’) 
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Fig.  Cll  -  ASYMPTOTIC  PEAK  PRESSURE  ISOBARS 


relative  amplitude 


Z  ,  ft 


Fig.  C 1 2  -  PEAK  DECAY  ALONG  Z  AXIS; INITIAL  VELOCITY 
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